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ABSTRACT 


The  Amplitude  Analysis  Programs  have  been  written  to  produce  an  average 
pulse  shape  over  a  reasonable  time  interval  for  up  to  200  range  stations 
behind  the  tracked  target  leading  edge  and  a  frequency  distribution  of  cross 
section  at  each  station  over  the  time  period.  The  final  output  is  presented 
as  Stromberg-Carlson  4020  plots  of  average  pulse  and  distributions.  The 
input  is  obtained  by  digitizing  A -Scope  films  on  the  Group  21  High  Precision 
computer  controlled  film  reader,  filtering  this  data  with  a  ramp  filter  to 
locate  the  position  of  the  tracked  target  leading  edge,  and  then  writing  a  tape 
containing  only  that  data  behind  the  leading  edge. 

The  Amplitude  Frequency  Distribution  and  Statistical  Program  uses  this 
filtered  output  tape  and  computes  an  average  pulse  and  the  standard  deviation 
from  this  pulse  at  each  desired  range  station  over  a  given  time  interval  and 
also  finds  the  distribution  of  cross  section  at  each  range  station,  if  desired. 

A  tape  containing  this  information  is  then  used  as  input  to  the  Amplitude 
Plotting  Program  which  produces  the  Stromberg-Carlson  4020  plots. 

The  present  report  contains  a  summary  of  the  method  of  obtaining  the  input 
data,  the  mathematics  involved,  a  description  and  listing  of  the  two  programs, 
and  directions  for  running  the  programs. 
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I.  INTRODUCTION 


The  Amplitude  Analysis  Programs  have  been  designed  to  produce  an  average 
pulse  shape  over  any  reasonable  time  interval  for  up  to  200  range  stations  along  a 
filmed  trace.  The  frequency  distribution  of  the  radar  backscattering  cross  section  at 
each  range  station  over  the  selected  time  period  is  also  available.  The  input  data  is 
obtained  from  the  digitization  of  filmed  A-Scope  data,  obtained  from  the  Wallops 
Island  Trailblazer  experiments,  on  the  Group  21  High  Precision  computer  controlled 
film  reader.  This  data  is  subsequently  filtered  using  a  ramp  filter  to  obtain  positions 
behind  the  tracked  target  leading  edge  before  being  entered  into  the  present  program. 
The  outputs  are  in  the  form  of  plots  of  average  pulse  shape  and  of  frequency 
distributions  of  radar  cross  section  at  each  range  station. 

11.  FILM  TRACE  DIGITIZATION 

The  Amplitude  Frequency  Distribution  and  Statistical  Program  accepts  a  digital 
definition  of  points  along  a  filmed  trace.  The  Group  21  High  Precision  computer 
controlled  automatic  film  reader  is  used  to  digitize,  pulse-by-pulse,  sets  of 
simultaneous  traces  from  A-scope  films.  Figure  1  shows  a  sample  film  to  be  digitized. 
The  missing  trace  seen  in  the  second  channel  corresponds  to  the  moment  of  recording 
a  time  mark  (see  channel  3).  Each  trace  will  produce  approximately  200  points. 

The  Y  axis  of  the  film  reader  is  visually  aligned  with  the  range  direction  of  the  traces 
so  that  points  of  common  Y  as  read  by  the  reader  can  be  considered  to  be  points  of 
common  range  for  each  trace  in  the  set.  The  traces  consist  of  either  quadrature 
pairs  or  independent  amplitude  recordings. 

The  output  of  the  film  reader  program  is  first  processed  by  the  "Filter”  program 
which  serves  to  edit  the  data  and  to  perform  a  range  tracking  function  in  identifying 
the  position  of  the  leading  edge  of  the  tracked  target.  This  computer  tracking  by 
means  of  the  pattern  recognition  of  a  pulse  and  the  identification  of  the  leading  edge 
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compensates  for  jitter  present  in  the  automatic  range  tracker  servo  system  of  the 
radar.  The  program  must  establish  a  zero  deflection  reference  for  each  trace  based 
only  on  the  data  describing  that  individual  trace.  The  pulse  return  generally  fills 
only  a  small  fraction  of  the  length  of  a  trace.  Therefore,  the  mode  of  the  distribution 
of  the  amplitude  values  of  each  trace  is  separately  calculated.  This  mode  is  used  as 
the  zero  reference  from  which  the  amplitude  and  polarity  for  each  point  of  that 
trace  are  computed.  Thus  the  data  are  self-fiducialized. 

The  ’’Filter"  program  is  used  to  obtain  the  leading  edge  of  the  tracked  target . 

The  amplitude  is  first  computed  either  directly  from  amplitude  data  or  by  translating 
the  bipolar  quadrature  data  to  a  unipolar  amplitude  function  by  summing  the  squares 
of  the  deflection  amplitudes  of  the  common  range  points.  The  amplitudes  are  then 
smoothed  using  a  ramp  filter  in  order  to  emphasize  the  leading  edge  of  a  target  pulse 
and  to  reduce  the  effect  of  erroneous  readings  caused  by  dirt  and  scratches  on  the 
film.  The  filtered  function  is  then  sliced  with  an  experimentally  determined  amplitude 
level  in  order  to  locate  the  position  of  the  tracked  target  leading  edge.  The  program 
then  outputs  the  corrected  data  for  this  leading  edge  position  and  all  subsequent 
increments.  In  this  way  returns  from  identifiable  range  stations  with  respect  to  the 
leading  edge  of  the  tracked  target  may  be  correlated  on  a  pulse-by-pulse  basis. 

The  output  from  the  Filter  program  is  an  800  bit  per  inch  packed  binary  tape, 

(Fig.  2).  The  first  record  in  each  file  is  the  BCD  title  from  the  input  (film  reader) 
tape.  Each  succeeding  record  contains  the  data  for  one  pair  of  traces.  The  first 
word  in  each  record  is  a  time  code,  the  first  18  bits  giving  the  time  block  number  and 
the  last  18  giving  the  number  of  the  current  trace  set  within  that  time  block.  Since 
a  time  trace  occurs,  in  general,  once  every  48  traces  (0.  05  sec. )  the  time  code  will 

be  interpreted  as  0-0,  0-1,  0-2, .  0-47,  1-0,  1-1, .  1-47,  etc.  When  running 

data  it  is  usual  to  start  on  a  time  trace  i.  e.  a  trace  whose  second  code  word  is  0 
(for  example,  0-0,  1-0,  2-0,  3-0,  etc.).  The  second  word  of  the  record  contains  a 
count  of  the  remaining  words  in  the  record  (i.  e.  the  number  of  recorded  points  behind 
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the  pulse  leading  edge).  The  third  and  fourth  words  are  not  used  by  this  program. 

The  remaining  words  in  each  record  are  each  divided  into  two  half-words  (18  bits): 
the  first  half-word  contains  channel  one  information;  the  second  half-word  contains 
channel  two  information.  In  the  case  that  only  one  channel  was  filtered,  the  second 
half-word  in  each  case  would  be  zero.  The  Amplitude  Frequency  Distribution  and 
Statistical  Program  is  equipped  to  unpack  these  words  and  assemble  them  in  usable 
form. 

III.  RADAR  CALIBRATION 

The  radar  system  is  calibrated  immediately  before  and  after  each  test.  The  entire 
system  is  calibrated  by  injecting  calibrated  test  signal  amplitudes  into  the  radar 
receiver  and  comparing  these  locally  generated  calibration  signals  with  the  signal 
amplitudes  returned  from  a  precision  metallic  sphere  of  known  cross  section  located 
at  known  distances  from  the  radar.  The  precision  sphere  is  carried  aloft  by  a 
balloon  of  negligible  radar  cross  section. 

This  sphere  and  signal -generator  calibration  enables  the  use  of  a  simplified 
radar  equation: 

PR  =  K  a/R4  (1) 


where 


P  -  power  received  in  milliwatts 

R 

R  -  range  in  nautical  miles 

a  -  radar  cross  section  in  square  meters 

4 

K  -  radar  coefficient  in  (nautical  miles)  milliwatts/meter 


Enough  information  is  obtained  from  this  sphere  and  signal -generator  combination  to 
determine  the  radar  system  coefficient  K  relating  radar  backscattering  cross  section 
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to  slant  range  and  signal -generator  calibrations.  The  constant  K  can  then  be  used  to 
compute  cross  sections  for  other  targets  whose  ranges  and  received  signal  amplitudes 
are  known. 

The  calibrated  signal  is  also  recorded  on  the  film  at  specified  power  increments 
between  -50  and  -115  dbm.  These  calibration  pulses  are  read  on  the  High  Precision 
automatic  film  reader  and  the  peak  displacement  from  baseline  at  each  power  level  is 
recorded.  A  calibration  curve  is  obtained  showing  power  in  dbm  as  a  function  of 
displacement  from  baseline.  A  sample  calibration  curve  is  shown  in  Fig.  3.  A  table 
of  values  of  power  in  dbm  vs.  amplitude  displacement  in  film  reader  units  is  entered 
into  the  Amplitude  Frequency  Distribution  and  Statistical  Program  to  be  used  in 
calibrating  the  data. 

IV.  MATHEMATICS 

The  Amplitude  Frequency  Distribution  and  Statistical  Program  is  designed  to 

produce  an  average  cross  section  over  a  given  time  interval  as  a  function  of  range 

station  behind  the  leading  edge  of  the  tracked  target  and  a  frequency  distribution  of 

the  cross  section  for  each  individual  station. 

The  raw  amplitude  data  (either  entered  directly  as  amplitude  from  the  input  tape 

or  computed  as  the  square  root  of  the  sum  of  the  squares  of  the  quadrature  data)  is 

2 

first  changed  to  cross  section  in  db/m  by  the  formula: 


a  (db)  =  40  log  R  -  K  (db)  +  DBM 


(-> 


where 


a  (db) 


radar  backscattering  cross  section  in  db 

2 

relative  to  one  square  meter,  (db/m  ). 


R 


radar  slant  range  in  nautical  miles. 
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K  (db) 


radar  coefficient  in  db 


DBM  -  received  signal  power  as  determined  from 

the  scaled  calibration. 

Then  the  cross  section  in  square  meters  is  found  by: 

/  \  <r(db)/10  /QX 

a  (square  meters)  =  10  (3) 

The  statistical  analysis  is  performed  on  the  cross  sections  in  square  meters  to 
give  an  average  pulse  and  the  standard  deviation  from  that  pulse  at  each  station.  The 
average  pulse  is  computed  by  the  equation 

N 

—  i  ~  1  °i  (Sfluare  meters) 

a  (square  meters)  = - — -  (4) 


where  N  is  number  of  points  over  which  averaging  is  performed.  The  standard 
deviation  is  then  calculated  as 


s.  d. 


\  1*1  <A  -2 

V - N - a 


(5) 


Any  signal  found  to  lie  above  the  calibration  curve  is  omitted  from  the  statistical 
analysis  since  its  true  power  value  can  not  be  determined;  any  signal  lying  below  the 
calibration  is  considered  as  noise  and  is  arbitrarily  set  equal  to  the  value  of  the 
minimum  discernible  signal,  (MDS). 

After  the  average  and  standard  deviations  have  been  computed  at  each  station  on  the 
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cross  section  in  square  meters,  they  are  converted  back  to  db/m  for  plotting. 
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A  frequency  distribution  may  also  be  taken  for  each  range  station  and  plotted  if 

desired.  This  distribution  is  found  by  simply  selecting  boxes  of  suitable  width  for 
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cross  section  in  db/m“  and  counting  the  number  of  points  at  each  station  which  fall 
into  each  box. 

V.  DESCRIPTION  OF  THE  AMPLITUDE  FREQUENCY  DISTRIBUTION  AND  STATISTICAL 
PROGRAM 

A.  General 

The  Amplitude  Frequency  Distribution  and  Statistical  Program  is  written  in 

Fortran  II  for  the  IBM  7094  computer  and  is  presented  in  Appendix  A.  The  program 

uses  as  input  the  packed  binary  tape  containing  the  digitized  traces  produced  by  the 

"Filter"  program.  The  major  output  is  a  BCD  800  bit  per  inch  tape  containing  a 

description  of  the  average  pulse  cross  sections,  the  average  pulse  ±  one  standard 

deviation,  and  a  frequency  distribution  of  cross  section  for  each  individual  station, 
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all  in  db/m~.  This  tape  is  used  as  input  to  the  Amplitude  Plotting  Program 

(see  Section  VI.).  Numerical  values  for  the  average  cross  section  and  standard 
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deviation  in  square  meters  and  in  db/m~  at  each  station  are  also  recorded  on  the 
monitor  tape. 

The  program  is  designed  to  handle  up  to  200  range  stations  per  trace  and  can  take 

a  distribution  over  up  to  100  distribution  boxes.  The  program  is  designed  to  make 

maximum  use  of  the  limited  core  storage  area  remaining.  A  generalized  flow  chart 

of  the  program  appears  in  Fig.  4.  The  current  record  is  first  read  from  the  input 

tape  and  converted  to  amplitude  in  film  reader  units  by  the  proper  subroutine  depending 

on  whether  it  is  desired  to  process  channel  one  alone,  channel  two  alone,  or  both 

channels  as  quadrature  data.  A  calibration  is  applied  to  convert  the  amplitude  in 

film  reader  units  to  cross  section  in  square  meters.  A  running  summation  of  cross 

section  and  of  cross  section  squared  at  each  station  is  computed  up  to  the  current 

/  2 

record.  The  cross  section  is  then  converted  back  to  db/m  and  counted  in  the  proper 
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distribution  box  for  that  value  of  cross  section  at  that  station.  The  next  record  is 
then  read  in  and  the  procedure  repeated  until  all  the  traces  to  be  processed  in  the 
current  set  have  been  read. 

When  all  the  traces  to  be  processed  in  the  present  set  have  been  read  into  the 

computer,  the  sum  of  the  cross  sections  and  the  sum  of  the  cross  sections  squared  for 

each  station  have  automatically  been  computed.  At  this  time  it  remains  only  to 

compute  average  cross  section  at  each  station  (Eq.  4)  in  square  meters  and  standard 

deviation  at  each  station  (Eq.  5).  The  average  and  the  average  ±  one  standard 

/  2 

deviation  are  then  computed  and  converted  to  db/m  ;  a  form  more  suitable  for  plotting. 
Pertinent  data  is  written  on  the  monitor  tape  and  the  averages  and  averages  ±  one 
standard  deviation  at  each  station  are  written  on  the  plotting  tape.  The  frequency 
distribution  for  each  station  is  also  written  on  the  plotting  tape  at  this  time. 

After  the  tapes  have  been  written  for  the  current  set  of  data,  the  program  will 
cycle  back  and  begin  to  process  the  next  set  as  above.  If  this  next  set  is  to  be 
independent  of  the  previous  set,  the  counters  are  cleared  and  the  program  re-run  with 
all  new  data.  If  it  should  be  desired  to  take  the  statistics  over  a  longer  time  interval, 
the  counters  are  not  cleared  at  this  time  and  the  next  set  of  data  is  simply  added  to  the 
previous  set.  The  total  number  of  traces  to  be  averaged  is  limited  only  by  possible 
overflow  of  the  summation  counters. 

B.  Subroutines 

The  Amplitude  Frequency  Distribution  Program  contains  7  Fortran  II  subroutines. 
Subroutine  AMPDA1  will  unpack  and  compute  amplitude  in  film  reader  units  using 
channel  one  data  only.  Subroutine  AMPDA2  will  unpack  and  compute  amplitude  using 
channel  two  data  only.  Subroutine  QUAD  will  compute  amplitude  using  both  channels 
as  quadrature  data  from  the  formula 

V2 

,22  22 

A  =  (A  sin  0  +  A  cos  0) 
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where  channel  one  is  assumed  to  contain  A  sin  0  and  channel  two  contains  A  cos  0 
(or  their  equivalent  in  polarized  amplitude  data). 

Subroutine  DBCOM  calibrates  the  data  i.  e.  substitutes  cross  section  in  square 
meters  for  the  given  amplitude  in  film  reader  units.  If  the  given  value  lies  above 
the  calibration  curve,  it  is  omitted  from  the  analysis  since  no  exact  measure  of  its 
power  can  be  found.  If  it  lies  below  the  curve,  the  power  is  arbitrarily  set  to  MDS 
for  that  point.  A  count  of  the  points  lying  above  and  below  the  curve  at  each  station 
is  kept  and  printed  out  on  the  monitor  tape.  A  linear  interpolation  is  used  to  find 
power  corresponding  to  amplitudes  lying  between  given  points  in  the  table.  The  cross 
section  in  square  meters  is  then  computed  by  Eqs.  (2)  and  (3).  Note:  since  for  ease 
in  punching,  the  -DBM's  of  the  calibration  are  read  in  as  positive  numbers,  Eq.  (2) 
must  be  rewritten  as 


cr  (db)  =  40  log  R  -  K  (db)  -  (-DBM) 

Subroutine  ADSTAT  simply  computes  the  sums  of  the  cross  sections  and  cross 
sections  squared  up  to  the  current  pulse  for  each  station  ignoring  any  point  lying  above 
the  calibration  curve. 

2 

Subroutine  FREQCT  changes  cross  section  back  to  db/m  and  adds  the  proper 
additive  factor  to  obtain  a  positive  number.  The  program  checks  to  make  sure  that 
the  value  lies  within  the  limits  of  the  boxes  (i.  e.  between  STBLOK  and  ENBLOK). 

If  the  value  lies  outside  these  limits,  a  special  code  counter  is  set.  If  the  value  lies 
within  the  limits,  the  count  in  the  proper  box  is  incremented  by  1  so  that  when  all 
traces  of  the  set  have  been  processed,  each  box  will  contain  the  frequency  of  values 
lying  within  that  box. 

Subroutine  STAT  computes  the  average  and  standard  deviation  at  each  station  over 
the  desired  time  interval.  The  number  of  points  to  use  is  first  determined  by 
subtracting  out  the  number  of  values  lying  above  the  calibration  curve  at  this  station 
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from  the  total  number  of  traces  used  in  this  set.  The  program  is  ended  here  if  any 
overflow  has  occurred  in  the  counters.  Averages  and  standard  deviations  at  each 
range  station  are  calculated  by  Eqs.  (4)  and  (5). 

The  program  also  uses  several  FAP  subroutines. 

Subroutine  SKIP  will  skip  NSKIP  records  within  one  file  on  the  input  tape  (A8) 
either  forward  if  NSKIP  is  positive  or  backward  if  NSKIP  is  negative. 

Subroutine  READA  will  read  the  next  binary  record  from  the  A8  tape  and  will 
unpack  the  time  code  and  the  number  of  data  points  in  the  record. 

Subroutine  UNPAK  is  used  to  unpack  the  data  words  which  consist  essentially  of 
two  18  bit  words  packed  as  one  36  bit  word.  The  first  18  bits  are  stored  as  a 
channel  one  reading  into  IXl  and  the  last  18  bits  are  stored  as  a  channel  two  reading 
into  1X2. 

C.  Outputs 

The  following  computed  outputs  are  written  on  the  monitor  tape  (A3)  for  each  range 
station. 

1.  Range  station  number. 

2.  Average  cross  section  in  square  meters. 

3.  Standard  deviation  in  square  meters. 

/  2 

4.  Average  cross  section  in  db/m  . 

.  2 

5.  Average  cross  section  4-  1  standard  deviation  in  db/m  . 

/  2 

6.  Average  cross  section  -1  standard  deviation  in  db/m  . 

7.  Number  of  points  used  in  analysis. 

8.  Number  of  points  lying  above  calibration  curve. 

9.  Number  of  points  lying  below  calibration  curve. 

The  number  of  points  lying  outside  the  limits  of  the  distribution  boxes  and  the  value  of 
the  lowest  cross  section  in  square  meters,  corresponding  to  MDS  are  also  recorded 
on  the  A3  tape. 
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The  other  output  tape  (B7)  is  written  at  800  bits/inch  BCD  and  contains  the  average 

cross  section  at  each  range  station  and  the  average  cross  section  ±  1  standard 
/  2 

deviation  all  in  db/m  .  If  desired  it  will  also  contain  the  frequency  distribution  of 
/  2 

cross  section  in  db/m  at  each  station.  The  tape  contains  the  following  information: 

I.  Frequency  distributions  consisting  of  one  group  of  data  for  each  station  as 
follows: 


Record  1  and  2 

Titles  to  print  on  plots. 

Record  3 

MDS  for  this  time  interval. 

Record  4 

2 

Average  cross  section  in  db/m  for  this  station. 

Records  5  -  end 

Value  of  cross  section  in  middle  of  distribution 
box  and  frequency  of  points  in  that  box  for  each 
of  the  EN  boxes. 

Format  (1  H  -  12  F  10.  3).  Each  record  thus 
contains  values  for  6  boxes. 

Data  for  the  next  station  follows  directly  in  the  same  format  until  all  stations 
have  been  recorded. 

II.  Average  Pulse  Cross  Section 

Each  time  interval  will  have  one  set  of  records  giving  the  computed  statistical 


data  as  follows: 

Record  1  and  2 

Titles  to  print  on  plots. 

Records  3— .  . . . 

Each  record  will  contain  data  for  1  station  in 

the  form: 

Station  number. 

Average  in  square  meters . 

Standard  Deviation  in  square  meters. 

Average  in  db/m^ 
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Average  +  1  standard  deviation  in  db/m  . 

2 

Average  -  1  standard  deviation  in  db/m  . 

Number  of  points  used  in  analysis. 

Number  of  points  above  calibration  curve. 

Number  of  points  below  calibration  curve. 

Format  (1H  -  15,  5F12.  5,  15,  214) 

After  sets  I  and  II  have  been  written  for  the  current  time  interval,  the  same  data 
for  the  next  time  interval  will  immediately  follow  on  the  tape.  If  no  frequency 
distribution  was  taken,  only  II  (the  average  pulse  cross  section  data)  will  appear  on 
the  tape.  Figure  5  shows  a  sample  diagram  of  the  arrangement  of  data  on  the 
B7  tape. 

VI.  THE  AMPLITUDE  PLOTTING  PROGRAM 

The  Amplitude  Plotting  Program  is  designed  to  produce  plots  of  the  frequency 

distributions  and  average  pulses  obtained  from  the  Amplitude  Frequency  Distribution 

and  Statistical  Program.  The  Plotting  Program  has  been  written  in  the  Fortran  IV 

language  and  is  run  under  the  IBSYS  monitor.  The  plots  are  produced  on  the 

Stromberg-Carlson  4020  plotter  and  are  available  either  as  hard  copy  or  film.  A 

flow  chart  of  the  program  is  shown  in  Fig.  6;  a  program  listing  is  given  in  Appendix  B. 

Figures  7a  and  7b  are  sample  frequency  distribution  plots.  It  will  be  seen  that 

the  data  from  two  stations  are  plotted  in  each  frame  and  that  the  titles  at  the  top  of  the 

frame  identify  the  data.  MDS  level  for  the  time  interval  is  shown  on  each  plot  and  the* 

average  signal  for  the  given  range  station  is  indicated.  The  distributions  are  plotted 

showing  number  of  cases  in  the  box  (or  frequency)  as  ordinate  as  a  function  of  cross 
2 

section  in  db/m  at  the  middle  of  the  box,  as  abscissa. 

Figure  7a  demonstrates  the  case  where  the  cross  section  distribution 
ranges  over  a  wide  set  of  values  corresponding  to  a  large  standard  deviation 
whereas  7b  demonstrates  the  case  where  the  distribution  is  much  more  compact 
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corresponding  to  a  small  standard  deviation.  The  peak  at  MDS  in  each  case  is 

influenced  by  the  arbitrary  cut-off  point  at  the  lower  end  of  the  calibration  curve 

below  which  all  data  is  considered  to  be  noise  and  set  arbitrarily  to  MDS  level. 

Figure  8  is  a  typical  plot  of  the  average  pulse  cross  section  taken  over  a  0.  1  sec. 

/  2 

interval.  The  value  of  cross  section  in  db/m  is  plotted  as  ordinate  against  range 

station  number  as  abscissa.  The  circles  indicate  the  position  of  the  average  cross 

section  taken  over  this  time  interval  at  the  indicated  station.  The  bars  indicate  the 

position  of  average  ±  one  standard  deviation.  The  standard  deviation  is  not  symmetric 

2 

about  the  average  when  plotted  on  the  db/m  scale  since  the  original  statistics  were 

performed  on  cross  section  in  square  meters  and  were  merely  changed  back  to 
2 

db/m  for  ease  in  plotting.  Thus,  although  one  standard  deviation  is  symmetric 
about  the  average  in  the  original  calculation,  it  is  not  symmetric  on  the  logarithmic 
scale  of  the  plots.  The  bottom  of  the  standard  deviation  is  arbitrarily  cut-off  at  MDS 
if  it  extends  below  this  level. 

VII.  OPERATING  PROCEDURES  FOR  THE  AMPLITUDE  FREQUENCY  DISTRIBUTION 
AND  STATISTICAL  PROGRAM 

The  Amplitude  Frequency  Distribution  and  Statistical  Program  is  run  under  the 
FORTRAN  monitor.  The  input  tape  is  the  800  bit/inch  binary  output  tape  from 
"FILTER"  and  is  mounted  on  A8.  If  a  tape  is  to  be  made  for  plotting  frequency 
distributions  or  average  pulses,  a  blank  tape  must  be  mounted  on  B7  and  set  to 
800  bit/inch.  This  tape  will  be  kept  and  used  as  input  to  the  Amplitude  Plotting  Program. 
The  program  and  control  data  input  is  from  A2  and  the  monitor  tape  is  A3.  Two  versions 
of  the  Amplitude  Frequency  Distribution  and  Statistical  Program  have  been  written. 
Version  I  is  the  program  described  above.  Version  II  is  a  simplified  version  of  the 
program  which  may  be  used  whenever  all  the  variables  in  the  program  remain 
constant  for  all  sets  of  data  and  where  it  is  desired  to  process  only  sequentially 
in  time,  i.  e.  whenever  it  is  desired  to  run  one  time  interval  and  then  run  the  next  set 
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by  repeating  the  last  time  trace  and  processing  the  next  sequential  time  interval 
immediately.  The  counters  are  cleared  at  this  time  so  that  each  set  must  be 
independent  of  the  previous  set.  The  actual  functioning  of  the  program  is  identical 
for  both  versions;  Version  II  is  usually  used  for  production  work  to  save  time  in 
punching  data  cards. 

The  control  cards  for  the  two  versions  of  the  program  are  made  out  as  follows: 


Control  Cards:  Version  I 

I.  The  first  card  following  the  "*DATA"  card  will  contain  6  fixed  point  numbers 
[Format  (615)] 

Name  Description  Columns 

NSET  The  number  of  sets  of  data  to  follow.  1-5 

Each  set  will  give  a  complete  analysis  of  any 
time  interval  desired.  Up  to  100  sets  of  data 
may  be  run  at  one  time:  however, 
each  set  must  be  in  the  same  file  on  the  input 
tape. 


IFB7  A  code  telling  if  averages  and  deviations  are  6-10 

to  be  written  on  B7  tape 

=  1  if  averages  and  deviations  are  to  be  written 
=  0  if  they  are  not  to  be  written 

NOFILE  Number  of  initial  files  on  the  A8  input  tape  to  11-15 

skip  before  beginning  processing — NOFILE  =  0 
for  first  file;  NOFILE  =  1  for  second  file,  etc. 
Usually  only  one  file  is  on  each  tape  so  that 
NOFILE  =  0. 
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Columns 


Name 

IFPLOT 


IXIN 


IFCAL 


Description 

Code  to  tell  if  frequency  distribution  is  to  be  16-20 
computed  and  written  on  B7 
=  1  if  distribution  to  be  written 
=  0  if  distribution  not  to  be  written 

Granularity  of  original  film  reader  data  21-25 

word — usually  IXIN  =  4  —  the  film 

reader  could  digitize  data  to  4000  counts 

per  inch  but  because  of  resolution  considerations 

only  1000  counts  per  inch  were  usually  utilized. 

Number  of  points  in  calibration  table  26-30 

=  0  if  no  calibration  to  be  used. 

=  number  of  points  in  table  if  calibration 
to  be  performed  —  0  <  IFCAL  <  50 
if  calibration  desired. 


The  second  data  card  will  be  a  title  card  punched  as  follows: 


Col.  1  must  contain  a  blank. 

Col.  2-5  contains  the  experiment  identification 
number  (ex.  A223). 


Col.  6-7  must  be  blank. 

Col.  8-72  any  title  desired  to  follow  the  time 
on  the  output  tape  title  —  usually 
starting  off  as  (in  Col.  8):  f,sec.  data  — ". 


The  third  data  card  will  contain  9  constants  to  be  used  throughout  all  sets 
[Format  (10F7.  5)]  Decimal  points  will  be  punched. 
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Name 

Description 

Columns 

XMULT 

Multiplicative  factor  for  channel  two. 

1-7 

XADD 

Additive  factor  for  channel  two. 

8-14 

YMULT 

Multiplicative  factor  for  channel  one. 

15-21 

YADD 

Additive  factor  for  channel  one. 

22-28 

AKLOG 

4  2 

Radar  constant  K,  in  db  (N.  M.  mw/m  ) 

29-35 

STBLOK 

Minimum  value  for  distribution  boxes 

—  STBLOK  >0. 

36-42 

ENBLOK 

Maximum  value  for  distribution  boxes. 

43-49 

EN 

Number  of  distribution  boxes.  (EN  <  100) 

50-56 

CONTPL 

Additive  constant  to  make  all  cross  section 

57-63 

values  in  db/m  positive  for  the  distribution 
determination. 


The  next  set  of  cards  will  have  ranges  punched  on  them  (one  for  each  set  to  be  analyzed). 
10  to  a  card — fFormat  (10F7.  0).  Decimal  points  will  be  punched.  The  ranges  will 
be  obtained  from  a  graph  of  range  vs.  time  or  from  the  output  of  the  Willmann  - 
Wallops  program.  Ranges  will  be  in  kilofeet. 

The  next  group  of  cards  will  contain  the  calibration,  if  one  is  to  be  used.  They 
will  be  punched  10  numbers  per  card — fFormat  ( 10F7 .  5)] .  Starting  with  the  first 
card  of  the  group  they  would  be  punched 


amplitude  1 

1-7 

-dbm  1 

8-14 

amplitude  2 

15-21 

-dbm  2 

22-28 

amplitude  3 

29-35 

-dbm  3 

36-42 

15 


amplitude  4 

43-49 

-dbm  4 

50-56 

amplitude  5 

57-63 

-dbm  5 

64-70 

Decimal  points  must  be  punched.  The  first  amplitude  recorded  will  be  the  lowest 
(highest  -dbm)  and  the  last  will  be  the  highest  amplitude  (lowest  -dbm). 

II.  Following  these  control  cards  which  are  fixed  for  the  entire  run,  will  follow 
NSET  sets  of  data  cards  of  2  cards  each  as  follows: 


Card  1  of  the  set  will  contain  8  fixed  point  numbers  —  Format  (1415). 


Name 

N 


M_ 

NSKIP 


IFCLFR 


Description  Columns 

Number  of  pulses  to  process  (as  many  as  1  -5 

desired). 

Last  station  per  pulse  to  process  (1  <  M  <200)  6-10 

Number  of  initial  records  to  skip  on  the  11-15 

A8  input  tape  (forward  or  back).  To  start 
at  the  beginning  of  the  file  NSKIP  =  0. 

Code  telling  whether  to  clear  the  frequency  16-20 
counters  and  sums  at  start  of  current  set 
=  1  if  desired  to  clear  counters 
=  0  if  desired  not  to  clear  counters  (in  general , 
the  counters  will  be  cleared  at  the  start  of  a 
set — however,  it  may  be  desired  to  continue 
averaging  over  a  longer  time  interval  in 
which  case  they  would  not  be  cleared.  For 
example,  to  compute  the  series  390.  0  -  390.  1, 
and  then  390.  0  -  390.  2  the  counters  would  not 
be  cleared  for  the  second  set). 
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Name 

NPT1 

I  NT 


Description  Columns 

Number  of  first  station  to  process.  21-25 

Interval  between  stations  to  process.  26-30 

i.  e.  INT  =  1  if  every  station  is  to  be 
processed;  INT  =  2  for  every  other 
station,  etc. 


IFSTAT 


IFQUAD 


The  second  card  of  the 
Punch  decimal  point. 


Code  telling  if  averages  and  standard 
deviations  to  be  computed. 

=  0  if  they  are  not  to  be  computed. 

=  1  if  they  are  to  be  computed. 

Code  telling  what  type  input  is  to  be  used 
=  -  1  if  only  channel  one  to  be  processed. 

=  + 1  if  only  channel  two  to  be  processed. 

=  0  if  both  channels  to  be  processed  as 
quadrature  data. 

set  will  contain  2  numbers  [Format  (2F10.  5)] 


31-35 


36-40 


Name  Description  Columns 

TMI  Initial  time  in  seconds  for  set.  1-10 

DTM  Interval  in  time  in  seconds  over  which  11-20 

Distribution  taken  (usually  DTM  =  0.  1). 

Version  II  has  been  written  to  make  the  program  as  simple  as  possible  to 
run  on  a  production  basis.  This  version  can  be  used  whenever  the  constants  do 
not  change  within  a  run  and  the  data  is  to  be  run  in  sequential  time  order. 
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Version  II  requires  the  following  data  cards. 


Card  1 


:  Format  ( J 415) 


Name 

NSET 


IFB7 


IFQUAD 


NOFILE 


IFPLOT 


Description 

Number  of  sets  of  data  to  process. 

Each  will  give  a  complete  analysis  of 
the  time  interval  desired.  Up  to  100  sets 
of  data  may  be  run  at  one  time  but 
must  be  sequential  in  time  as  described 
above. 

Code  telling  if  averages  and  deviations 
to  be  written  on  B7  tape. 

=  1  if  averages  and  deviations  to  be  written 
=  0  if  they  are  not  to  be  written 

Code  telling  what  type  input  is  to  be  used 
=  -1  if  only  channel  one  to  be  processed. 

=  +1  if  only  channel  two  to  be  processed. 

=  0  if  both  channels  to  be  processed  as 
quadrature  data. 

Number  of  initial  files  on  the  A8  input  tape  to 
skip  before  beginning  processing — NOFILE 
for  first  file;  NOFILE  =  1  for  second  file, 
etc. 

Code  to  tell  if  frequency  distribution  to  be 
computed  and  written  on  B7. 

=  1  if  distribution  to  be  written. 

=  0  if  distribution  not  to  be  written. 


Columns 

1-5 


6-10 


11-15 


16-20 

0 


21-25 
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Name 

IXIN 


IFCAL 


N_ 

M_ 

NSKIP 

IFCLER 

NPT1 

INT 


Description 

Granularity  of  original  film  reader  data 
word — usually  IXIN  =  4 — the  film  reader 
could  digitize  data  to  4000  counts  per  inch 
but  because  of  resolution  considerations  only 
1000  counts  per  inch  were  usually  utilized. 

Number  of  points  in  calibration  table. 

=  0  if  no  calibration  to  be  used. 

=  number  of  points  in  table  if  calibration  to 
be  performed  — 0  <  IFCAL  <  50  if 
calibration  to  be  used. 

Number  of  pulses  to  process  (as  many  as 
desired). 

Last  station  to  process  per  pulse  (1  <M<200) 

Number  of  initial  records  to  skip  on  input 
tape  A 8 — can  only  skip  records  before  the 
first  set  for  this  version. 

Must  be  +1  in  this  version  as  counters  must 
be  cleared. 

Number  of  first  station  to  process. 

Interval  between  stations  to  process,  i.  e. 

INT  =  1  to  process  every  station;  INT  =  2  for 
every  other  station,  etc. 


Columns 

26-30 


31-35 


36-40 

41-45 

46-50 

51-55 

56-60 

61-65 
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Name 
IF  ST  AT 


Description 

Code  telling  if  averages  and  standard 
deviations  to  be  computed. 

=  0  if  not  to  be  computed. 

=  1  if  they  are  to  be  computed. 


Columns 

66-70 


Card  2  will  be  a  title  card  punched  as  follows: 


Col.  1  must  contain  a  blank. 

Col.  2-5  contain  the  experiment  identification 

number  (A223  for  example). 

Col.  6-7  are  blank. 

Col.  8-72  any  title  desired  to  follow  the  time  on  the 

output  —  usually  starting  off  in  Col.  8 
with  the  words  "sec.  data  —  ". 


Card  3:  [Format  (2F10.  5)]  Punch  decimal  points. 


initial  time  of  run,  in  seconds. 

interval  in  time  in  seconds  over  which  to  take 

distribution  (usually  =  0.  1). 


1-10 

11-20 


Card  4:  [Format  (10F7.  5)]  Punch  decimal  points. 

Description 


Mame 

XMULT 

XADD 

YMULT 

YAPP 

AKLOG 


Columns 


Multiplication  factor  for  channel  two. 

Addition  factor  for  channel  two. 

Multiplication  factor  for  channel  one. 

Addition  factor  for  channel  one. 

4  2 

Radar  constant  K  in  db  (N.  M.  mw/m  ) 


1-7 

8-14 

15-21 

22-28 

29-35 
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Name  Description  Columns 

STBLOK  Value  of  minimum  distribution  box  36-42 

(0  <  STBLOK). 

ENBLOK  Maximum  distribution  box.  43-49 

EN  Number  of  boxes  (EN  <  100).  50-56 

CONTPL  Additive  constant  to  make  all  values  of  57-63 


cross  section  in  db/m  positive  for  the 
distribution  determination. 

The  next  group  of  cards  will  have  ranges  punched  on  them  (one  for  each  set  computed). 
They  will  be  punched  10  to  a  card  in  Format  (10F7.  0).  Decimal  points  will  be  punched. 
The  ranges  will  be  obtained  from  a  graph  of  range  vs.  time  for  each  0.  1  sec.  interval 
or  from  the  output  of  the  Willmann  -  Wallops  Program.  Ranges  will  be  in  kilofeet. 

The  last  set  of  cards  will  contain  the  calibration  used.  They  will  be  punched 
as  follows  10  numbers  per  card  in  Format  (10F7.  5)  starting  with  the  first  card. 


amplitude 

i 

1-' 

7 

-dbm  1 

8- 

15 

amplitude 

2 

15 

-21 

-dbm  2 

22 

-28 

amplitude 

3 

29 

-35 

-dbm  3 

36 

-42 

amplitude 

4 

43 

-49 

-dbm  4 

50 

-56 

amplitude 

5 

57 

-63 

-dbm  5 

64 

-70 
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The  first  amplitude  recorded  will  be  the  lowest  (highest  -dbm)  and  the  last  will  be  the 
highest  amplitude  (lowest  -dbm). 

VIII. OPERATING  PROCEDURE  FOR  AMPLITUDE  PLOTTING  PROGRAM 


The  Amplitude  Plotting  Program  will  plot  the  frequency  distributions  and  the 
average  pulses  on  the  Stromberg-Carlson  4020.  It  uses  as  input  the  B7  tape  obtained 
from  the  Amplitude  Frequency  Distribution  and  Statistical  Program  (which  will  be 
mounted  on  B7)  and  produces  as  output  a  4020  tape  on  B6.  The  program  is  run 
under  1BSYS. 

Input  consists  of  two  cards  for  each  set  of  data  with  common  constants  (usually 
the  entire  run  may  be  made  without  changing  constants). 

Card  1:  Format  (515) 


Name 

MM 

NBLOK 


M 


IF  FREQ 


Description  Columns 

Number  of  time  intervals  to  process.  1-5 

Number  of  boxes  over  which  distribution  was  6-10 
ta  ken . 

Number  of  stations  per  time  interval  to  plot.  11-15 
—  if  not  every  station  is  processed,  M 
here  is  only  the  number  of  actual  stations 
run  and  not  the  last  station  run. 

=  0  if  frequency  distribution  to  be  plotted  16-20 

=  1  if  distribution  not  to  be  plotted  and  does 
not  exist  on  the  input  tape. 

=  -1  if  no  distribution  is  to  be  plotted  but  the 
distribution  has  been  written  on  the  input  tape. 
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Card  2:  Format  (14F5.0)  Punch  decimal  points. 


Name 

XMINFR 

XMAXFR 

DXFR 

YMAXFR 

DYFR 

XMINA 


XMAXA 

DXA 

YMINA 

YMAXA 

DYA 


Description 

Minimum  limit  of  distribution  boxes 

2 

(minimum  abscissa)  in  db/m~. 

Maximum  limit  of  distribution  boxes 

2 

(maximum  abscissa)  in  db/m  . 

Increment  for  tic  marks  on  x  axis  of 

2 

distribution  in  db/m  . 


Columns 

1-5 


6-10 


11-15 


Maximum  frequency  to  plot  (for  0.  1  sec.  16-20 

intervals  YMAXFR  =  100)  (maximum  ordinate). 


Increment  for  tic  marks  on  y  axis  of  21-25 

distribution. 

Minimum  x  scale  value,  minimum  station  26-30 
number,  for  average  pulse  plot — usually 
XMINA  =  0. 

Maximum  x  scale  value,  maximum  station  31-35 
number,  for  average  pulse  plot. 

Increment  in  x  axis  tic  marks  for  average  36-40 
pulse  plot. 


Minimum  y  scale  value  in  db/m"-  for  average  41-45 
pulse  plot. 

,  2 

Maximum  y  scale  value  in  db/m  for  average  46-50 
pulse  plot. 

/  2 

Increment  in  y  axis  tic  marks  in  db/m  for  51-55 
average  pulse  plot. 


23 


The  program  will  produce  the  frequency  distributions  of  two  stations  per  frame. 
Hence,  if  the  number  of  stations  (M)  is  even,  the  frequency  distributions  will  occupy 
M/2  frames  per  time  interval.  If  M  is  odd,  they  will  occupy  M/2  +  1  frames  (since 
the  extra  station  must  be  plotted  somewhere).  Also,  there  is  a  title  frame  and  an 
average  pulse  frame  for  each  time  interval.  Thus,  the  total  number  of  frames 
required  should  be  (if  frequency  distributions  are  plotted): 

For  M  even  -  number  of  frames  =  MM  (M/2  +  2) 

or 

For  M  odd  -  number  of  frames  =  MM  (M/2  +  3) 

If  only  average  plots  are  to  be  plotted,  there  will  be  only  MM  frames. 
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Y  AXIS 


3-21-6790 


CHANNEL  3 - 

CHANNEL  2 - 1 

CHANNEL  1 - 1 


X  AXIS 


Fig.  i.  Sample  A-Scope  Film. 
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3-21-6791 


Somple  Input  Tope 
Format 


BCD  Title 

Troce  Set  1 

Troce  Set  2 

Troce  Set  3 

• 

• 

• 

Troce  Set  N 

EOF 

BCD  Title 

Troce  Set  1 

Trace  Set  2 

Troce  Set  3 

• 

• 

• 

Troce  Set  M 

EOF 

2(o) 


Somple  Doto  Record 
Format  for  A  Typical 


Set  of 

Traces 

IB  Bits^ 

I  B  BJts 

Time 

Block 

1  Troce 

J  In  Set 

K 

1 

|  Unused 

1 

1 

Unused 

1 

1 

Unused 

1 

Pt  1 
Chon  1 

Pt  1 

J  Chon  2 

Pt  2 

Chon  1 

1  Pt  2 

J  Chon  2 

Pt  3 
Chon  1 

J  Pt  3 
[  Chon  2 

• 

• 

• 

Pt  K 
Chon  1 

Pt  K 

[  Chon  2 

END  OF 
RECORD  GAP 


} 


2(b) 


Word  1 

Word  2 
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Fig.  2.  Sample  Input  Tape  Format- 


3-21 -6795 


Fig.  3.  Typical  Calibration  Curve  Relating  Power  in  -DBM 
to  Amplitude  in  Film  Reader  Units. 
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3-21-6792 


Fig.  4.  Flow  Chart  of  Amplitude  Frequency  Distribution 
and  Statistical  Program. 
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Fig.  5.  Sample  Output  Tape  Format. 
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3-21-6793 


Fig.  6.  Flow  Chart  of  Amplitude  Plotting  Program  . 
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105 

4223  3*9. 9- 400,0  SEC  DMA  FILTERED  M  26$  6  5fA  BEFORE  LE 
AMPLITUDE  FROH  *  31 N  PHI  AND  4  COS  PHl 
HlHl HUH  DISCERNIBLE  3I6HAL - -1 10. GG  CB.H 

STATION  no.  23  AVERAGE  SIGNAL-  -15-303  £DB  ABOVE  1  U .  ME  TER I 


100. 


Fig.  7(a).  Sample  Frequency  Distribution  with  Large  Standard 
Deviation. 
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Fig.  7(b). 


Sample  Frequency  Distribution  with  Small  Standard 
Deviation. 
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A223  399. *-400.0  SEC  OAT*  F1LTERE0  AT  268  6  STA  BEFORE  tE 
AMPL1TU0E  FROM  A  SIN  PHI  AND  A  COS  PHI 
AVERAGE  PU.SE  CROSS  SECTION 
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Fig.  8.  Sample  Plot  of  Average  Pulse  Cross  Section  Showing  One 
Standard  Deviation  of  Either  Side  of  Average. 
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WILSON ♦FLORENCE  AMPLITUDE  PROGRAM  VERSION  1  11/5/65 

LIST 

LABEL 

SYMBOL  TAbLE 

DIMENSION  STATEMENTS 

DIMENSION  A  <  50 ) ♦ AMP< 2  20 )  ♦ AVER ( 200 )  * Db  <  50  )  > I C T < 200 ♦ 1O4 ) > I  TEMP  <  220 ) 
If  IX  1220)  ♦  S  ANUM  (200)  *SD(200)  ♦SUMAMP(200)  tSUMSO<200) 

DIMENSION  AVDBM { 200 )  ,AVDBP<200) ♦AVDBI200)  ♦  X  1< 200  )  ♦  X4 ( 200  ) 

DIMENSION  TITLE< 12) *RNG< 100) 

COMMON  STATEMENTS 

COMMON  I T 1 ♦ I T2  .NOP ♦ EN ♦ENBLOK ♦ I  NT ♦  M ♦ NB LOK ♦ NP T 1 ♦STBLOK ♦ W I DTH ♦ I  X  IN , EM 
1 ♦ ANUM ♦ CON TPL  fNCAL  fCONLOG 

COMMON  XMULTtXADD, YMULT ♦ YADD  ♦ AMPLOW ♦ DBLOW ♦ AMPH I ,DBHI 
COMMON  IX ♦ I  TEMP ♦ AMP ♦ AVER ♦ I CT ♦ SD ♦ SUMAMP ♦ SANUM ♦ Db ♦ A ♦ SUMSQ 
FORMAT  STATEMENTS 
FORMAT { 12A6 ) 

FORMAT ( 6H  T I  ME = I  5 ♦ 3X ♦ I  5 ) 

FORMAT { 1H1 ) 

FORMAT { A6 ♦ F  5 • 1 ♦ 1 H-F 5 . 1 ♦ 1 1 A6 ) 

FORMAT { 39H  AMPLITUDE  FROM  A  SIN  PHI  AND  A  COS  PHI) 

FORMAT ( 1415 ) 

FORMAT ( 72H 
1  ) 

FORMAT { 3H  N = I  4 ♦ 4X ♦ 2HM = I  4 ♦ 4X ♦ 6HNSK I P= I  4 ♦ 4X ♦ 7H I F C LER = I  4 ♦ 4X ♦ 6H I F C AL = I 
14t4Xt7HlFOUAD=I4) 

FORMAT { 7F 1 0 • 5 ) 

FORMAT { lHU4XtlHJ9X.8HAVER  MSQ12X.6HSD  MS07X ♦ 5 HA VDbM8X ♦ 4HDBH I  8  X  >  4HD 
lBL05Xt 1HN2X.2HHI 2X.2HLO) 

FORMAT { 1 2H  STATION  NO . I  4 ♦ 6 X ♦ 1 5HA VE R AGE  S I GN AL  =  F8 . 3 ♦ 23H  < Db  ABOVE  1 
1  SQ.  METER ) ) 

FORMAT (lXI5t 1P2E18«8*0P3F12*5»I5«2I4) 

FORMAT { 39H  AMPLITUDE  FROM  CHANNEL  1  ONLY  ) 

FORMAT { 8H  S T BLOK  = F 8 . 2 ♦ 7 X ♦ 7 HE NBLOK= F 8 . 2 ♦ 4X ♦ 2HN = F 5 . 0 ♦ 5X ♦ 6HW I DTH = F 8 . 3 


C 

C 

c 

c 

c 

c 

c 

c 


1 ) 

FORMAT ( 1H  12F10.3) 

FORMAT ( 39H  AMPLITUDE  FROM  CHANNEL  2  ONLY  ) 

FORMAT ( 3H  X=F10*4t2H  +F10.4.6H  Y=F10.4.2H  +F10.4) 

FORMAT { 28H  MINIMUM  DISCERNIBLE  S I GNAL  =  F  8 • 2  *4H  DBM) 

FORMAT { 1H  I5»5F12«5*I5»2I4) 

FORMAT { 1  OF  7 . 5 ) 

FORMAT { 7H0MAX I M= I  5 ) 

FORMAT ( 7H  R ANGE  =  F 10 *4 ♦ 5X ♦ 2HK  =  F 1 0 . 4/ 19H  CROSS  SECT  I ON  =  DBM+F10 • 5/22H 
1  LOWEST  CROSS  SEC T I  ON  =  F 1 2 . 8  ♦  1 4H  SOUARE  METERS) 

FORMAT { 1 H  6F 1 5 • 5 ) 

FORMAT { 20H0CAL I  BRAT  I  ON  FOLLOWS) 

FORMAT { lH04XtlHJ6Xt6HAVMESQ7X ♦ 5HSDMSQ6X ♦ 5HA VDBM8X ♦ 4HDBH I 8X ♦AHDBLOS 
1 X  >  1HN2X  ♦2HHI2Xt2HLO) 

CONVERSION  FROM  KILOFEET  TO  NAUTICAL  MILES 
AALOG=40#*LOG1 OF (6.0761033 ) 

NSET  =  NUMBER  OF  TIME  INTERVALS  TO  PROCESS 
I FB7= 1  IF  AVERAGE  AND  DEVIATION  TO  BE  WRITTEN  ON  B7 

I FB7=0  IF  NO  AVERAGE  AND  DEVIATION  TO  BE  WRITTEN  ON  B7 

N0FILE=NUMBER  OF  INITIAL  FILES  TO  SKIP  ON  INPUT  TAPE 
I F PLOT  =  1  IF  FREOUENCY  DISTRIBUTION  TO  BE  WRITTEN  ON  B7 

I FP LOT  =  0  IF  NO  FREOUENCY  DISTRIBUTION  TO  BE  WRITTEN  ON  B7 

IXIN= GRANULARITY  OF  ORIGINAL  READING 
I FC AL=NUMBER  OF  CALIBRATION  POINTS  IN  CURVE 
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C  I FCAL=0  IF  NO  CALIBRATION  TO  BE  PERFORMED 

READ  INPUT  TAPE  2 » 1 02  ♦  NSE T , I FB7 , NOF I LE , I FPLOT , I  X  I N , I F CAL 
READ  INPUT  TAPE  2, 197, TITLE 

C  XMULT=MULT IPL ICAT ION  FACTOR  FOR  CHANNEL  2 

C  XADD=ADD I T I VE  CONSTANT  FOR  CHANNEL  2 

C  YMULT  =  MULT IPLICAT ION  FACTOR  FOR  CHANNEL  1 

C  YADD  =  ADD  I T I VE  CONSTANT  FOR  CHANNEL  1 

C  AKLOG=RADAR  K  FACTOR 

C  STBLOK=MINIMUM  END  OF  FREQUENCY  BLOCKS 

C  ENBLOK=MAXIMUM  END  OF  FREQUENCY  BLOCKS 

C  EN=NUMBER  OF  BLOCKS 

C  CONTPL=ADDI T IVE  CONSTANT  TO  MAKE  ALL  VALUES  POSITIVE 

C  FOR  DISTRIBUTION 

READ  INPUT  TAPE  2 , 1 1 6 , XMUL T  , X ADD , YMUL T , Y ADD , AKLOG , S TB LOK , ENBLOK , 
1 EN , CON  TPL 

C  RNG=RANGE  OF  THIS  SET  IN  KFT 

READ  INPUT  TAPE  2,116, (RNG(J)  »J=1*NSET) 

CALL  DNSITYI 1,8,2) 

READ  IN  CALIBRATION  TABLE  IF  CALIBRATION  TO  BE  DONE 
A=AMPL I TUDES  FOR  CALIBRATION  TABLE 
DB=POWER  IN  -DBM  FOR  CALIBRATION  TABLE 
I  F ( IFCAL)92»93»92 
92  NCAL=IFCAL 

READ  INPUT  TAPE  2 , 1 1 6  ,  ( A ( J J  )  , DB ( J J )  , J J  =  1  , NC AL 
AMPLOW=SMALLEST  AMPLITUDE  POSSIBLE 

DBLOW  =  POWER  CORRESPONDING  TO  MINIMUM  AMPLITUDE  READING 
AMPH I =H I GHEST  AMPLITUDE  POSSIBLE 

DBH I =POWER  CORRESPONDING  TO  MAXIMUM  AMPLITUDE  READING 
AMPLOW=A ( 1 ) 

DBLOW=DB ( 1 ) 

AMPH I = A ( NCAL  ) 

DBHI =DB ( NCAL ) 

NEXT  A3  PAGE 
RATION  TABLES  ON  A3 


93 

4 

5 
76 


75 

6 


SKIP 

TO  TOP  OF 

WRITE 

OUTPUT 

TAPE 

3*99 

WRIT 

E  OUT 

CALI 

WRITE 

OUTPUT 

TAPE 

3*12 

WRITE 

OUTPUT 

TAPE 

3*12 

SET 

UP  B7 

TAPE 

I  F ( I FB7  )  4 , 5 , 
REWIND17 

4 

CALL 

DNSITYI 

2*7**:  ! 

1 

I , ( A ( J J ) , DB ( J J ) , J J= 1 , NCAL 
IF  IT  IS  TO  BE  USED 


IF ( I FPLOT 176,75,76 
CALL  DNSITY(2»7»2) 

REWIND  17 

SHUT  OFF  INDICATORS  AT  START  TO  CHECK  LATER 
IF  ACCUMULATOR  OVERFLOW  6,6 
IF  DIVIDE  CHECK  7,7 

FIND  PROPER  FILE  ON  A8  INPUT  TAPE 

REWIND  8 

CALL  AHEAD  (8, NOFILE, 0) 

READ  BCD  TITLE  FROM  INPUT  TAPE 
READ  INPUT  TAPE  8,103 

PERFORM  LOOP  FOR  EACH  SET  OF  DATA 
DO  50  I SET= 1 »NSET 

N=NUMBER  OF  TRACES  TO  PROCESS 
M=NUMBER  OF  LAST  STATION  TO  PROCESS 
NSK I P=NUMBER  OF  INITIAL  RECORDS  TO  SKIP 

IFCLER=1  IF  COUNTERS  TO  BE  CLEARED  FOR  THIS  SET  OF  DATA 
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C  I FCLER=0  IF  COUNTERS  NOT  TO  BE  CLEARED  AT  THIS  TIME 

C  NPT 1=NUMBER  OF  FIRST  STATION  TO  PROCESS 

C  I N T  = I NTERVAL  BETWEEN  STATIONS  TO  PROCESS 

C  I FST AT  =  1  IF  STATISTICAL  ANALYSIS  TO  BE  PERFORMED 

C  I  FSTAT  =  0  IF  NO  STATISTICAL  ANALYSIS  TO  BE  PERFORMED 

C  I FQUAD=0  FOR  QUADRATURE  DATA 

C  IFQUAD=NEGATIVE  IF  CHANNEL  1  ONLY  TO  BE  USED 

C  IFQUAD=POSITIVE  IF  CHANNEL  2  ONLY  TO  BE  USED 

READ  INPUT  TAPE  2 , 1 02  . N  . M  ,  NSK I P  ,  I FC L ER , NP T 1 . 1  NT . I F ST A T  *  I  FQUAD 
TM1=START ING  TIME  FOR  THIS  SET 

DTM  =  T I  ME  INTERVAL  OVER  WHICH  SET  TO  BE  PROCESSED 
READ  INPUT  TAPE  2 , 1 05 , TM 1 . DT M 


303  TM2=TM1+DTM 


WRITE 

OUT 

TITLES  AND  CODE  WORDS  ON  A3  TAPE  FOR  REFERENCE 

WRITE 

OUTPUT 

TAPE 

3*99 

WRITE 

OUTPUT 

TAPE 

3*100*TITLE(1)*TM1*TM2*(TITLE(J)*J=2*12) 

I  F ( I FQUAD ) 400 

»  40 1 

*402 

400 

WRITE 

OUTPUT 

TAPE 

3*109 

GO  TO 

403 

402 

WRITE 

OUTPUT 

TAPE 

3*112 

GO  TO 

403 

401 

WRITE 

OUTPUT 

TAPE 

3*101 

403 

WRITE 

OUTPUT 

TAPE 

3*103 

WRITE 

OUTPUT 

TAPE 

3  *  1 04  *  N *M  *  NSK I P  *  I FC L ER *  I F C AL *  I FGU AD 

FIND 

PROPER  POINT  ON  INPUT  TAPE  TO  BEGIN  PROCESSING 

CALL  SKIP(NSKIP) 

CLEAR  COUNTERS  IF  DESIRED 
ALSO  SET  SUMS  SO  FAR  TO  ZERO 
I F ( I  FC  LER  )  9  *  1 0  *  9 
9  DO  26  1=1*200 
DO  11  J  =  1  *  1 04 
ICT ( I , J ) =0 
11  CONTINUE 

SUMAMP ( I ) =0. 

SUMSQ ( I ) =0. 

26  CONTINUE 
ANUM=0. 

C  GET  CURRENT  RANGE 

10  RANGE=RNG ( ISET) 

C  COMPUTE  WIDTH  OF  COUNTING  BLOCK 

W I D  TH= ( ENBLOK-STBLOK ) /EN 

WRITE  OUTPUT  TAPE  3 » 1 1 0 , STBLOK , ENBLOK , EN » W I DTH 
WRITE  OUTPUT  TAPE  3  *  1 1 3 » XMU L T » XADD , YMUL T , Y ADD 
C  C0NL0G=40  LOG  R  -  K 

RANL0G=L0G1 OF (RANGE ) 

CCLOG=AKLOG+AALOG 

CONLOG=40«*RANLOG-CCLOG 

C  CSLOW=MINIMUM  DISCERNIBLE  SIGNAL 

CSLOW=-DB( 1 l+CONLOG 
CSLOW=10.** ( CSL0W/10. ) 

WRITE  OUTPUT  TAPE  3* 122  .RANGE  , AKLOG » CONLOG » CSLOW 
C  NN = AN UM= NUMBER  OF  POINTS  OVER  WHICH  STATISTICS  TO  BE  TAKEN 

ANUM=ANUM+FLOATF ( N ) 

NN= ANUM 
NBLOK=  EN 
DO  51  JK  =  1  * N 

C  READ  NEXT  RECORD  FROM  A8 

CALL  READA 
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C  CHECK  TIME  CODE  IF  FIRST  RECORD  OF  SET 

IF (JK-1 ) 300*300,301 

300  WRITE  OUTPUT  TAPE  3.98.IT1.IT2 

301  NOP=NOP 

C  REVERSE  WORDS  INTO  FORTRAN  ORDER  IN  MATRIX 

DO  13  J=1 *NOP 
JJ=22 1-J 
ITEMPI J)=IX( JJ> 

13  CONTINUE 

C  COMPUTE  AMPLITUDE  BY  PROPER  METHOD  DEPENDING  ON  TYPE  DATA 

I F ( IFQUAD) 14,15,16 

14  CALL  AMPDA1 
GO  TO  17 

16  CALL  AMPDA2 
GO  TO  17 

15  CALL  QUAD 

C  CALIBRATE  DATA 

17  I F ( IFCAL)90,18,90 

90  CALL  DBCOM 

C  PERFORM  ADDITION  FOR  THIS  RECORD 

18  CALL  ADSTAT 

C  DO  FREQUENCY  DISTRIBUTION  IF  DESIRED 

I F ( I FPLOT ) 9 1 ,51 ,91 

91  CALL  FREQCT 

C  END  OF  LOOP - GO  READ  NEXT  RECORD  IF  ANY  REMAINS  IN  SET 

51  CONTINUE 

END  OF  LOOP - READY  TO  DO  STATISTICAL  ANALYSIS  OVER  PERIOD 

COMPUTE  AVERAGE  AND  STANDARD  DEVIATION  EACH  STATION 
IF ( I FSTAT ) 19,20,19 

19  CALL  STAT 
WRITE  OUTPUT  TAPE  3,106 

COMPUTE  AVERAGE, AVERAGE  +  ONE  S.D.  AND  AVERAGE  -  ONE  S.D. 
CHANGE  BACK  TO  DB  OVER  1  SQUARE  METER 
DO  52  J=NPT1,M,INT 
NNN=  SANUM ( J ) 

AVDBl J)=10.*LOG10F(AVER( J) ) 

AVDBPI J)=AVER I JJ+SDIJ) 

AVDBPI J > =10.*LOG10F < AVDBP < J)  ) 

AVDBMI J ) =AVER ( J )-SD( J ) 

IF  BELOW  MDS  SET  TO  VALUE  AT  MDS 
I F  t AVDBMI JJ-CSLOWI2 20,2 20, 221 

220  AVDBMI J > =CSLOW 

221  AVDBMI J ) =10. *LOG10F ( AVDBMI J ) ) 

WRITE  VALUES  ON  MONITOR  TAPE  A3 

ALSO  SHOW  IF  ANY  POINTS  ABOVE  OR  BELOW  CALIBRATION  CURVE 

222  WRITE  OUTPUT  TAPE  3 , 1 08 , J , AVER  I J ) , SD I J )  , A VDB ( J )  , AVDBP I J  )  , A VDBM I J )  , 
1NNN, ICT I J, 103 ) , ICT I J, 104) 

52  CONTINUE 

C  COMPUTE  MAXIMUM  NUMBER  IN  ANY  BLOCK  AND  PRINT  VALUE 

20  I F I  I FPLOT )  21,22,21 

21  MAXI M=  0 

DO  70  J=NPT 1 ,M ♦ I  NT 
DO  70  J J J  =  1 ♦ NBLOK 
IF! ICT (J,JJJ)-MAXIM)70,70,71 
71  MAXIM=ICT (J,JJJ) 

70  CONTINUE 

WRITE  OUTPUT  TAPE  3, 120, MAXIM 
C  WRITE  4  LINES  OF  TITLES  ON  PLOTTING  TAPE 
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WRITE  OUTPUT  TAPE  17 , 1 00  ,  T I T LE (1)  . TM1  ,  TM2  ,  < T I TLE < J )  ,  J  =  2 , 1 2  ) 

I  F (  IFQUAD)404»405»406 

404  WRITE  OUTPUT  TAPE  17.109 
GO  TO  407 

406  WRITE  OUTPUT  TAPE  17,112 
GO  TO  407 

405  WRITE  OUTPUT  TAPE  17,101 

407  ALA=-DB(1) 

WRITE  OUTPUT  TAPE  17, 114, ALA 
XI  ( 1  )  =STBL0K+WIDTH/2.  -CONTPL 
NBLM1=NBL0K-1 

C  STORE  MID  POINT  EACH  BLOCK  IN  XI 

DO  82  JKL=1,NBLM1 
L=JKL+1 

X 1 ( L ) =X 1 ( JKL1+WIDTH 

82  CONTINUE 
DO  72  J  =  NP  T 1 »M  » I  NT 

WRITE  OUTPUT  TAPE17,107,J  ,AVD6(J) 

DO  83  JKL=1 ,NBLOK 
X4(JKL)=ICT ( J , JKL ) 

83  CONTINUE 

WRITE  MID  POINT  OF  BLOCK  AND  FREQUENCY  IN  THAT  BLOCK 
ON  PLOTTING  TAPE 

202  WRITE  OUTPUT  TAPE  17 . 1 1 1  ,  < X 1 < J J >  ,  X4 < J J )  ,  J J= 1 , NbLOK ) 

72  CONTINUE 

C  WRITE  VALUES  FOR  AVERAGE  PULSE  ON  PLOTTING  TAPE 

22  I  F ( I FB7 ) 24 , 50 , 24 

24  WRITE  OUTPUT  TAPE17.100  ,  T I TL E ( 1 )  ,  TM 1  ,  TM2  ,  ( T I TL E ( J )  ,  J  =  2 , 1 2  ) 

IF ( I FQUAD)408»409,410 

408  WRITE  OUTPUT  TAPE  17,109 
GO  TO  411 

410  WRITE  OUTPUT  TAPE  17,112 
GO  TO  411 

409  WRITE  OUTPUT  TAPE17.101 

411  DO  53  J  =NP T 1  ,M , I N T 
NNN=  SANUM ( J ) 

WRITE  OUTPUT  TAPE  1 7 , 1 1 5  ,  J  ,  A VER ( J )  ,  SD ( J )  , A VDB ( J )  , AVDBP ( J ) , AVDBM ( J ) 
1  ,NNN »ICT ( J , 1 0  3 ) »ICT ( J , 104 ) 

53  CONTINUE 

C  GO  DO  NEXT  SET 

50  CONTINUE 

C  REWIND  AND  UNLOAD  TAPES  AT  END 

I F ( I FB7 128,29,28 

28  END  FILE  17 
CALL  RUN (2,7) 

GO  TO  208 

29  I  F ( I FPLOT  )  30,208 ,30 

30  END  FILE  17 
CALL  RUN (2,7) 

208  CALL  EXIT 
END 
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DIMENSION  STATEMENTS 

DIMENSION  A (50)  ,  AMP (220)  ,AVER(200)  *DB(50)  *  I C T ( 200  *  1 04  )  *  I  TEMP (220) 
*  I  X ( 220  )  *  S ANUM (200)  ,SD(200)  ,SUMAMP(200)  ,SUMSQ(200) 

DIMENSION  AVDBM ( 200 ) , AVDBP(200) ,AVDB(200) ,X1(200)*X4(200) 

DIMENSION  T  I  T L E ( 12 )  ,  RNG  ( 100) 

COMMON  STATEMENTS 

COMMON  I  Tin  T2  ,  NOP,EN,ENBLOK, I  NT  ,M,NBLOK,NPT 1 ,STBLOK* WIDTH, IX  IN, EM 
,  ANUM, CON TPL,NCAL,CONLOG 

COMMON  XMULT , X ADD, Y MULT , Y  ADD  ,  AMPLOW , DBLOW , AMPH I  ,DBHI 
COMMON  1 X, I  TEMP  *  AMP , AVER , I CT  *  SD  *  SUM AMP , SANUM, DB  *  A  *  SUMSQ 
FORMAT  STATEMENTS 
FORMAT ( 12A6 ) 

FORMAT ( 6H  T I ME= I  5 , 3X , I  5 ) 

FORMAT ( 1 H 1 ) 

FORMAT ( A6,F5* 1 , 1H-F5. 1 , 11A6) 

FORMAT ( 39H  AMPLITUDE  FROM  A  SIN  PHI  AND  A  COS  PHI) 

FORMAT (1615) 

FORMAT ( 7  2  H 

) 

FORMAT (3H  N= I  4 , 4X  ,  2HM= I  4 , 4X , 6HNSK I P= I  4 , 4X , 7H I FCLER= I  4 , 4X , 6H I FCAL= I 
4,4X*7HIFQUAD=I4) 

FORMAT ( 7F10.5 ) 

FORMAT ( 1H04X, 1HJ9X, 8HAVER  MSQ12X,6HSD  MSQ7X , 5HAVDBM8X , 4HDBH I  8 X  * 4HD 


BL05X, 1HN2X»2HHI 2X,2HLO) 

FORMAT ( 1 2  H  STATION  NO  .  I  4 , 6 X  ,  1 5HAVERAGE  S I GNAL  =  F  8 • 3 , 2  3H  ( DB  ABOVE  1 
SQ.  METER) ) 

FORMAT (lXI5tlP2E18.8,0P3F12.5*I5,2I4) 

FORMAT ( 3 9H  AMPLITUDE  FROM  CHANNEL  1  ONLY  ) 

FORMAT ( 8H  STBLOK=F8.2 , 7X , 7HENBLOK=F8 • 2 , 4X , 2HN*F5 • 0 , 5X , 6HW I DTH=F8 • 3 
) 

FORMAT ( 1H  12F10.3) 

FORMAT (39H  AMPLITUDE  FROM  CHANNEL  2  ONLY  ) 

FORMAT ( 3H  X=F1Q.4,2H  +F1G.4,6H  Y=F1G.4,2H  +F1Q.4) 

FORMAT ( 2 8H  MINIMUM  DISCERNIBLE  S I GNAL=  F 8.2* 4H  DBM) 

FORMAT (1H  I5,5F12.5*I5,2I4) 

FORMAT ( 10F7.5 ) 

FORMAT ( 7H OM AX  I M= I  5 ) 

FORMAT ( 7H  RANG E=F 1 0 . 4 , 5X , 2HK=F 1 0 • 4 / 1 9H  CROSS  SECT  I ON  =  DBM  +  F 10 . 5/22H 
LOWEST  CROSS  SECT  I ON  =  F 1 2 • 8 , 1 4H  SQUARE  METERS) 

FORMAT ( 1 H  6F15.5) 

FORMAT (20H0CALIBRAT ION  FOLLOWS) 

FORMAT ( 1H04X,1HJ6X  *6HAVMESQ7X , 5HSDMSQ6X , 5HAVDBM8X ,4HDBHI8X,4HDBL05 
X ♦ 1HN2X  *  2HH I 2X , 2HLO ) 

CONVERSION  FROM  KILOFEET  TO  NAUTICAL  MILES 


A ALOG =4 0.*LOG1 OF ( 6.0761033 ) 

NSET  =  NUMBER  OF  TIME  INTERVALS  TO  PROCESS 
I FB7= 1  IF  AVERAGE  AND  DEVIATION  TO  BE  WRITTEN 


ON 


I FB 7  =  0  IF  NO 
I F QU AD  =  Q  FOR 
I FQUAD=NEGAT I VE 
IFQUAD=POSIT ive 


AVERAGE  AND  DEVIATION 
QUADRATURE  DATA 

IF  CHANNEL  1  ONLY  TO 
IF  CHANNEL  2  ONLY  TO 


TO  BE  WRITTEN 


B7 

ON 


B7 


BE 

BE 


USED 

USED 
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nnnn  non  n  nnnnnnnnnn  n n 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 


NOF I LE=NUMBER  OF  INITIAL  FILLS  TO  SKIP  ON  INPUT  TAPt 
I FPLOT  =  1  IF  FREQUENCY  DISTRIBUTION  TO  BE  WRITTEN  ON  B7 
I FPLOT  =  0  IF  NO  FREQUENCY  DISTRIBUTION  TO  BE  WRITTEN  ON  B7 
IXIN= GRANULARITY  OF  ORIGINAL  READING 
I FCAL=NUMBER  OF  CALIBRATION  POINTS  IN  CURVE 
I FC AL=0  IF  NO  CALIBRATION  TO  BE  PERFORMED 
N=NUMBER  OF  TRACES  TO  PROCESS 
M=NUMBER  OF  LAST  STATION  TO  PROCESS 
NSK I P=NUMBER  OF  INITIAL  RECORDS  TO  SKIP 

I FC LER  =  1  IF  COUNTERS  TO  BE  CLEARED  FOR  THIS  SET  OF  DATA 
I FCLER=0  IF  COUNTERS  NOT  TO  BE  CLEARED  AT  THIS  TIME 
NPT 1 =NUMBER  OF  FIRST  STATION  TO  PROCESS 
I NT= I  NT  ERVAL  BETWEEN  STATIONS  TO  PROCESS 
I  FSTAT  =  1  IF  STATISTICAL  ANALYSIS  TO  BE  PERFORMED 
I FST AT  =  U  IF  NO  STATISTICAL  ANALYSIS  TO  BE  PERFORMED 
READ  INPUT  TAPE  2  *  1 02 * NSET *  I FB7 *  I FQUAD * NOF I LE *  I FPLOT *  I  X  I N *  I FCA L 
1 *N*M*NSKIP* IFCLER*NPT1 *  I  NT *  I FST AT 
READ  INPUT  TAPE  2,197*TITLE 

TM1  =  ST  ART  I NG  TIME  FOR  THIS  SET 

DTM=  T I  ME  INTERVAL  OVER  WHICH  SET  TO  BE  PROCESSED 
READ  INPUT  TAPE  2 , 105  * TM1 *DTM 

XMULT=MULTIPLICATI0N  FACTOR  FOR  CHANNEL  2 
X ADD= ADD I T I VE  CONSTANT  FOR  CHANNEL  2 
YMULT=MULTIPLICAT ION  FACTOR  FOR  CHANNEL  1 
Y ADD  =  ADD  I T I VE  CONSTANT  FOR  CHANNEL  1 
AKLOG=RADAR  K  FACTOR 

STBLOK=MIN IMUM  END  OF  FREQUENCY  BLOCKS 
ENBLOK  =  MAX I  MUM  END  OF  FREQUENCY  BLOCKS 
EN  =  NUMBER  OF  BLOCKS 

CONTPL=ADDITI VE  CONSTANT  TO  MAKE  ALL  VALUES  POSITIVE 
FOR  DISTRIBUTION 

READ  INPUT  TAPE  2  *  1 16  * XMUL T * X ADD  * YMUL T * Y ADD  * A KLOG * ST B LOK * E NBLOK * 

1 EN  *  CON  T  PL 

RNG=RANGE  OF  THIS  SET  IN  KFT 
READ  INPUT  TAPE  2  *  1 16  * ( RNG ( J )  * J  =  1 *NSET ) 

CALL  DNSITY ( 1*8*2  ) 

READ  IN  CAL  I  BRAT  ION  TABLE  IF  CALIBRATION  TO  BE  DONE 
A  =  AMPL I TUDES  FOR  CALIBRATION  TABLE 
DB=POWER  IN  -DBM  FOR  CALIBRATION  TABLE 
I  F ( I FCAL ) 92  *  93  *  92 
92  NC AL= I FCA  L 


93 

4 


READ  INPUT  TAPE  2  *  1 16  * ( A ( J J ) * DB ( JJ )  * JJ= 1  * NCAL  ) 

AMPLOW= SMALLEST  AMPLITUDE  POSSIBLE 

DBLOW  =POWER  CORRESPONDING  TO  MINIMUM  AMPLITUDE  READING 
AMPH I =H I GHEST  AMPLITUDE  POSSIBLE 

DBH I =POWER  CORRESPONDING  TO  MAXIMUM  AMPLITUDE  READING 
AMPLOW  =  A ( 1 ) 

DBLOW=DB( 1 ) 

AMPH I =A ( NCAL ) 

DBHI=DB(NCAL) 

SKIP  TO  TOP  OF  NEXT  A3  PAGE 
WRITE  OUTPUT  TAPE  3*99 

WRITE  OUT  CALIBRATION  TABLES  ON  A3 
WRITE  OUTPUT  TAPE  3*124 

WRITE  OUTPUT  TAPE  3  *  1 2 3  * < A < J J ) , DB < J J )  ,  J J= 1  * NC AL  ) 

SET  UP  B7  TAPE  IF  IT  IS  TO  BE  USED 
I  F (  IFB7)4,5*4 
REWIND17 
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n> 


c 


c 


c 

c 

c 


c 

c 


c 


c 

c 


c 


c 


c 


CALL  DNS  I T Y ( 2  *  7  *  2  ) 

5  I  F  (  IFPL0TI76, 75*76 
76  CALL  DNSITY(2»7»2) 

REWIND  17 

SHUT  OFF  INDICATORS  AT  START  TO  CHECK  LATER 
75  IF  ACCUMULATOR  OVERFLOW  6.6 

6  IF  DIVIDE  CHECK  7,7 

FIND  PROPER  FILE  ON  A8  INPUT  TAPE 

7  REWIND  8 

CALL  AHEAD  (8, NOFILE, 0) 

READ  BCD  TITLE  FROM  INPUT  TAPE 
READ  INPUT  TAPE  8,103 

PERFORM  LOOP  FOR  EACH  SET  OF  DATA 
DO  50  ISET=1 »NSET 

BACKSPACE  ONE  RECORD  AFTER  FIRST  SET 
I F  t I  SET-1  1302,303,302 


OF  INTERVAL 


302 

NSKIP 

=  -l 

GET  T 

IME 

at  end 

303 

TM2=TM1+DTM 

WRITE 

OUT 

TITLES 

WRITE 

OUTPUT 

TAPE 

3,99 

WRITE 

OUTPUT 

TAPE 

3,100, 

I  F  t  IFQUAD1420 

*  42  1 

,422 

420 

WRITE 

OUTPUT 

TAPE 

3,109 

GO  TO 

423 

422 

WRITE 

OUTPUT 

TAPE 

3,112 

GO  TO 

423 

421 

WRITE 

OUTPUT 

TAPE 

3,101 

423 

WRITE 

OUTPUT 

TAPE 

3,103 

WRITE 

OUTPUT 

TAPE 

3,104, 

FIND 

PROP 

ER  POIN 

CALL  . 

SKIPtNSK 

IP) 

AND  CODE  WORDS  ON  A3  TAPE  FOR  REFERENCE 


1 )  ,  TM 1  ,  TM2  ,  (TITLE (J) »J  =  2»12) 


1  1 


CLEAR  COUNTERS 
ALSO  SET  SUMS  SO 
DO  26  1=1,200 
DO  11  J= 1 , 104 
ICT  t  I  ,  J  )  =0 
CONTINUE 


ON  INPUT  TAPE  TO  BEGIN  PROCESSING 


FAR  TO  ZERO 


SUMAMP  t I  ) =0. 

SUMSQt I ) =0. 

26  CONTINUE 
ANUM  =  0  • 

GET  CURRENT  RANGE 
10  RANGE  =  RNG  t ISET ) 

COMPUTE  WIDTH  OF  COUNTING  BLOCK 
W I DT  H=  t ENBLOK-STBLOK ) /EN 

WRITE  OUTPUT  TAPE  3 » 1 1 0 , STBLOK , ENBLOK , EN , W 1 DTH 
WRITE  OUTPUT  TAPE  3,113, X MULT  ,XADD,YMULT .YADD 
CONLOG=40  LOG  R  -  K 
RANLOG=LOG10F( RANGE) 

CCLOG=AKLOG+AALOG 

CONLOG=40.*RANLOG-CCLOG 

CSLOW=MINIMUM  DISCERNIBLE  SIGNAL 
CSLOW=-DB  t 1 ) +CONLOG 
CSLOW=10.**(CSLOW/10.  ) 

WRITE  OUTPUT  TAPE  3 , 1 22 , RANGE , AKLOG , CONLOG , CSLOW 

NN=ANUM=NUMBER  OF  POINTS  OVER  WHICH  STATISTICS  TO  BE  TAKEN 
ANUM=ANUM+FLOATF(N) 
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n  n  n  r»  n  n  n 


NN=ANUM 
nblok=en 
DO  51  JK  =  1  * N 

C  READ  NEXT  RECORD  FROM  A8 

CALL  READA 

C  CHECK  TIME  CODE  IF  FIRST  RECORD  OF  SET 

I F ( JK-1 ) 300.300*301 

300  WRITE  OUTPUT  TAPE  3.98.IT1.IT2 

C  REVERSE  WORDS  INTO  FORTRAN  ORDER  IN  MATRIX 

301  NOP=NOP 

DO  13  J  =  1  *  NOP 
J J=22 1 - J 
ITEMPI J  )  =  I  X ( JJ) 

13  CONTINUE 

C  COMPUTE  AMPLITUDE  BY  PROPER  METHOD  DEPENDING  ON  TYPE  DATA 

I  F (  IFOUAD) 14.15.16 

14  CALL  AMPDA1 
GO  TO  17 

16  CALL  AMPDA2 
GO  TO  17 

15  CALL  QUAD 

C  CALIBRATE  DATA 

17  IF ( IFCAL)90,18,90 

90  CALL  DBCOM 

C  PERFORM  ADDITION  FOR  THIS  RECORD 

18  CALL  ADSTAT 

C  DO  FREQUENCY  DISTRIBUTION  IF  DESIRED 

I F ( IFPLOT ) 9 1 .51 .91 

91  CALL  FREQCT 

C  END  OF  LOOP - GO  READ  NEXT  RECORD  IF  ANY  REMAINS  IN  SET 

51  CONTINUE 

END  OF  LOOP - READY  TO  DO  STATISTICAL  ANALYSIS  OVER  PERIOD 

COMPUTE  AVERAGE  AND  STANDARD  DEVIATION  EACH  STATION 
I  F (  IFSTAT ) 19.20.19 

19  CALL  STAT 
WRITE  OUTPUT  TAPE  3.106 

COMPUTE  AVERAGE  .AVERAGE  +  ONE  S.D.  AND  AVERAGE  -  ONE  S.D. 
CHANGE  BACK  TO  DB  OVER  1  SQUARE  METER 
DO  52  J=NPT1 *M, INT 
NNN=SAMUM( J ) 

AVDBl J)=10.*LOG1UF( AVER! J) ) 

AVDBP ( J ) =AVER ( J ) +SD( J ) 

AVDBP! J)=10.*LOG10F( AVDBP! J) ) 

AVDBM(J)=AVER(J)-SD(J) 

IF  BELOW  MDS  SET  TO  VALUE  AT  MDS 
IF ( AVDBM ( J ) -CSLOW 1220,220,221 

220  AVDBM ( J ) =CSLOW 

221  AVDBM! J ) =10. *LOG10F ( AVDBM! J ) ) 

WRITE  VALUES  ON  MONITOR  TAPE  A3 

ALSO  SHOW  IF  ANY  POINTS  ABOVE  OR  BELOW  CAL  I  BRAT  I  ON  .CURVE 

222  WRITE  OUTPUT  TAPE  3 » 108 » J » AVER ( J ), SD ( J ), A VDB ( J ), AVDBP ( J ), AVDBM ( J )  , 
1NNN,  ICT (J.1C3) *ICT(J*104) 

52  CONTINUE 

COMPUTE  MAXIMUM  NUMBER  IN  ANY  BLOCK  AND  PRINT  VALUE 

20  IF(IFPLOT)  21,22,21 

21  MAX  I M  =  0 
DO  70  J=NPT1 ,M, INT 
DO  70  J J J= 1 » NBLOK 
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n 


c 


I F ( ICT (J,JJJ> -MAXIM) 70,70,71 

71  MAXIM=ICT( J.JJJ) 

70  CONTINUE 

WRITE  OUTPUT  TAPE  3, 120, MAXIM 

WRITE  4  LINES  OF  TITLES  ON  PLOTTING  TAPE 
WRITE  OUTPUT  TAPE  1  7 , 1 00  ,  T  I  TL E ( 1  )  , TM 1 , TM2  ,  ( T I TLE ( J )  , J  =  2  ,  1 2  ) 

I F (  I FQUAD1400 ,401 ,402 
400  WRITE  OUTPUT  TAPE  17,109 
GO  TO  403 

402  WRITE  OUTPUT  TAPE  17,112 
GO  TO  403 

4 U 1  WRITE  OUTPUT  TAPE  17,101 

403  ALA=-Db(l) 

WRITE  OUTPUT  TAPE  17, 114, ALA 
XI ( 1 ) =STbLOK+WIDTH/2.  -CONTPL 
NBLM1=NBL0K-1 
C  STORE  MID  POINT  EACH  BLOCK  IN  XI 

DO  82  JKL  = 1 , NBLM1 
L  =  JKL+  1 

X  1  ( L ) =X 1 ( JKL ) +WIDTH 

82  CONTINUE 
DO  72  J  =NP T 1 , M , I N T 

WRITE  OUTPUT  TAPE17,107,J  ,AVDB(J) 

DO  83  JKL  = 1 »NBLOK 
X4( JKL)=ICT( J.JKL) 

83  CONTINUE 

WRITE  MID  POINT  OF  bLOCK  AND  FREQUENCY  IN  THAT  BLOCK 
ON  PLOTTING  TAPE 

202  WRITE  OUTPUT  TAPE  1  7  ,  1 1 1  ,  (  X 1  (  J  J  )  ,  X4  (  J  J  )  ,  J  J  =  1 ,  NBJ-OK  ) 

72  CONTINUE 

C  WRITE  VALUES  FOR  AVERAGE  PULSE  ON  PLOTTING  TAPE 

22  I  F (  IFB7)24»350»24 

24  WRITE  OUTPUT  TAPE17.100  ,  T I TL E ( 1 )  , T M 1 , T M2  ,  ( T I TL E ( J ) , J  =  2 , 1  2  ) 

IF ( IFQUAD)404,405,406 

404  WRITE  OUTPUT  TAPE  17,109 
GO  TO  407 

406  WRITE  OUTPUT  TAPE  17,112 
GO  TO  407 

405  WRITE  OUTPUT  TAPE17.101 

407  DO  53  J=NPT1,M,INT 
NNN=SANUM ( J ) 

WRITE  OUTPUT  TAPE  1 7 , 1 1 5  ,  J  ,  A VER ( J ) , SD ( J  )  , AVOB ( J )  , AVDBP ( J ) , AVDBM ( J ) 
1 ,NNN, ICT ( J, 103 ) » ICT( J.104) 

53  CONTINUE 

C  COMPUTE  STARTING  TIME  FOR  NEXT  INTERVAL 

350  TM1=TM1+DTM 

C  GO  DO  NEXT  SET  IF  ANY 

50  CONTINUE 

C  REWIND  AND  UNLOAD  TAPES  AT  END 

550  IF(IFB7)28,29,28 

28  END  FILE  17 
CALL  RUN ( 2 , 7  ) 

GO  TO  208 

29  I F ( IFPLOT 130,208,30 

30  END  FILE  17 
call  run (2,7) 

208  CALL  EXIT 
END 
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*  LIST 

*  LABEL 

*  SYMBOL  TABLE 

SUBROUTINE  AMPDA1 

AMP  1  AMPLITUDE  USING  CHANNEL  1  ONLY 
DIMENSION  STATEMENTS 

DIMENSION  A (50) .AMP (220) ♦ AVER ( 200 )  * DB ( 5 0  )  .  I C T ( 200 . 104 )  .ITEMPI 220) 
1* I  X ( 2  20)  *SANUM ( 200  )  .SDI200)  ,SUMAMP(2  00)  .SUMSQ <2 00) 

C  COMMON  STATEMENTS 

COMMON  ITl»IT2»NOP»EN,ENULOK»INT»M»NbLOK»NPTl»STBLOK»WIDTH»IXIN»EM 
1  *  ANUM  *  CONTPL *NCAL ,CONLOG 

COMMON  XMULT .XADD.YMULT » YADD  .AMPLOW  .DBLOW  .AMPH I  ,DBHI 
COMMON  IX  *  I  TEMP. AMP .AVER , I C T , SD , SUMAMP , SANUM , DB  .  A . SUMSQ 
C  UNPACK  AND  COMPUTE  AMPLITUDE  FOR  CHANNEL  1  ONLY 

DO  1  J=NPT 1 »M* I  NT 
CALL  UNPAKI  ITEMPI J)  .1X1. 1X2) 

AMP ( J ) = I  X 1 / I  X  I N 
1  AMP ( J) =AMP ( J ) *YMULT+YADD 
RETURN 
END 


*  LIST 

*  LABEL 

*  SYMBOL  TABLE 

SUBROUTINE  AMPDA2 

CAMP2  AMPLITUDE  USING  CHANNEL  2  ONLY 
C  DIMENSION  STATEMENTS 

DIMENSION  A (50) , AMP (220)  .AVER (200)  .DBI50) . ICT ( 200. 104 ) , I  TEMP (220) 
1.IXI220)  .SANUM (200)  .SDI200)  .SUMAMP (200)  .SUMSQ (200) 

C  COMMON  STATEMENTS 

COMMON  1T1.1T2,N0P.EN,ENBL0K,INT,M,NBL0K.NPT1.STUL0K.WIDTH.IX1N,LM 
1 , ANUM. CONTPL ,NCAL .CON LOG 

COMMON  XMULT .XADD.YMULT .YADD  . AMPLOW , DbLOW .AMPH 1 , DBH I 
COMMON  I  X . I  TEMP , AMP .AVER. ICT , SD . SUMAMP  , SANUM , DB  . A . SUMSQ 
C  UNPACK  AND  COMPUTE  AMPLITUDE  FOR  CHANNEL  2  ONLY 

DO  1  J=NPT1 »M» INT 
CALL  UNPAKI  ITEMPI J) .1X1. 1X2) 

AMP( J)=iX2/IXIN 
1  AMP( J)=AMP( J ) *XMULT+XADD 
RETURN 
END 


*  LIST 

*  LABEL 

*  SYMBOL  TABLE 

SUBROUTINE  QUAD 

CQUAD  AMPLITUDE  FROM  QUADRATURE  DATA  USING  BOTH  CHANNELS 
C  DIMENSION  STATEMENTS 

DIMENSION  A ( 50 )  ,AMP( 220),AVER(200).DB(50).ICT(200.104),  I  TEMP ( 220 ) 
1  .  1X1220)  .SANUM ( 200  )  ,SD( 200)  .SUMAMP ( 200 )  .SUMSQ! 200 ) 

C  COMMON  STATEMENTS 

COMMON  ITI. IT2.NOP.EN . ENBLOK. INT *M .NbLOK ,NPT 1 .STBLOK. WIDTH, IXIN.LM 
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1 ,ANUM,CONTPL,NCAL,CONLOG 

COMMON  XMULT.XADD.YMULT.YADD  , AMPLOW , DBLOW . AMPH I , OBH 1 
COMMON  IX  ,  I  TEMP , AMP .AVER , I  CT  » SD  , SUM AMP , SANUM , DB . A . SUMSO 
C  UNPACK  AND  COMPUTE  AMPLITUDE  FROM  BOTH  CHANNELS 

DO  1  J  =  NPT1 *M» I  NT 
CALL  'JNPAK  (  I  TEMP  (  J  )  ,  I  XI  .  1X2  ) 

3  AM 1=1X1  / I  X  I N 
AM2 -  I  X 2/ I  X  I N 
AM1=AM1*YMULT+YADD 
AM2=AM2*XMULT+XADD 
5  AMP( J) =SQRTF ( AM1*AM1+AM2*AM2 > 

1  CONTINUE 
RETURN 
END 
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n  n  n  n 


*  LIST 

*  LABEL 

*  SYMBOL  TABLE 

SUBROUTINE  DBCOM 

CDBCOM  ROUTINE  TO  CALIBRATE  THE  DATA 
C  DIMENSION  STATEMENTS 

DIMENSION  A (50)  ,AMP(220) , AVER < 200 ) » DB < 50 )  ,  I  CT ( 200  ♦  104 )  .  I  TEMP (220) 
1.IXI220)  .SANUM ( 200 )  .SD(200)  .SUMAMP(200)  .SUMSQ1200) 

C  COMMON  STATEMENTS 

COMMON  IT1 , IT2.N0P.EN.ENBL0K, I  NT .M.NBLOX.NPT 1 .STBLOK .WIDTH. IXIN.LM 
1  , ANUM.CONTPL  »NCAL .CON LOG 

COMMON  XMULT  .XADD.YMULT . YADD  . AMPLOW , DBLOW . AMPH  I  , DBH I 
COMMON  IX. I  TEMP. AMP. AVER. ICT . SD . SUM AMP , SANUM , DB . A . SUMSQ 
C  PERFORM  LOOP  FOR  EACH  STATION 


DO  1  J  =  NPT 1 *M* I  NT 

C  IF  VALUE  ABOVE  CURVE  COUNT  IT  AS  BAD  POINT 

IF(AMP(J)-AMPHI 12.2.3 

3  ICT(J,103)=ICT(J. 1031+1 

C  ALSO  SET  AMP  =  777 77  (CODE) 

AMP(J) =77777. 

GO  TO  1 

C  IF  BELOW  CURVE  SET  VALUE  EQUAL  TO  LOWEST  POWER 

2  I F ( AMP ( J ) -AMPLOW 14,5.5 

4  AMP ( J ) =DBLOW 

C  AND  COUNT  AS  POINT  BELOW  CURVE 

ICT( J,104)  =  ICT( J,  1041+1 
GO  TO  8 

5  DO  12  1=1, NC AL 

INTERPOLATE  TO  FIND  POWER  CORRESPONDING  TO  GIVEN  AMPLITUDE 
BY  LINEAR  INTERPOLATION  BETWEEN  TABULATED  POINTS 
I F ( AMP ( J ) -A ( I  )  ) 10,11  ,  12 
11  AMP ( J ) =DB ( I ) 

GO  TO  8 

10  AAB= ( AMP ( J ) —A ( I  — 1 )  ) /  (  A (  I  ) — A ( I  — 1 )  ) 


DDB= ( DB ( 1 ) -DB ( I -1 ) ) * AAB 
14  AMP ( J ) =DB ( 1-1 ) +DDB 
GO  TO  8 
12  CONTINUE 

DBM  READ  IN  AS  POSITIVE  NUMBER  THUS  MUST  TAKE  -DBM  READ 
CROSS  SEC  T I  ON ( DB ) = -D6M+40  LOG  R  -K 
8  AMP ( J ) =-AMP ( J ) +CONLOG 
C  CROSS  SECTION  IN  SQUARE  METERS 

AMP ( J ) =10.** ( AMP ( J ) / 10. ) 

1  CONTINUE 
RETURN 
END 
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n  n 


*  LIST 

*  LAbEL 


AVERAGES 


♦  I  TEMP  (220  1 


*  SYMBOL  TABLE 

SUBROUTINE  ADSTAT 

CADST  SUBROUTINE  TO  ADD  AND  STORE  FOR 
C  DIMENSION  STATEMENTS 

DIMENSION  A (50) ,AMP(220) , AVER ( 200 ) . Db ( 5 0 ) . I C T ( 200 . 104 1 
1  *  IX ( 220  ) * SANUM ( 200 )  .50(200)  .SUMAMP ( 200 )  .SUMSQt  200 ) 

C  COMMON  STATEMENTS 

COMMON  I T 1  *  I T2  *  NOP  .EN.ENBLOL , I N T ,M , NBLOK , NP T 1 , STBLOK . W I D TH , I  X  I N  ,  LM 
1 .ANUM.CONTPL.NCAL  *  CON LOG 

COMMON  XMULT  .XADD.YMULT *YADD  , AMPLOW » DbLOW . AMPH I » DBH I 
COMMON  IX. I  TEMP.  AMP  .AVER.  I C T  .  SD  .  SUMAMP  .  SANUM , Db  .  A . SUMSQ 
C  IF  ABOVE  CURVE  DO  NOT  USE  IN  STATISTICS 

DO  1  J  =  NPT 1 *M  *  I  NT 
I F ( AMP ( J ) -77777 • ) 10 . 1 » 1 

FIND  SUM  OF  CROSS  SECTION 
UP  TO  CURRENT  POINT 
10  SUMAMP ( J  )  =  SUMAMP ( J ) +AMP ( J ) 

SUMSQ ( J  )  =  SUMSG ( J)+AMP(J)*AMP(J) 

1  CONTINUE 
RETURN 
END 


AND  CROSS  SECTION  SQUARED 
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n  n  n  n  n  n  n 


*  LIST 

*  LABEL 

*  SYMBOL  TABLE 

CFREQ  FREQUENCY  DISTRIBUTION 

SUBROUTINE  FREQCT 

c  dimension  statements 

DIMENSION  A (50)  ,AMP(220) . AVER ( 200 ) »DB ( 50 ) * ICT ( 200 . 104 ) . I  TEMP (220) 
1.IXI220)  .SANUM ( 200 ) »SD(200) ,SUMAMP<200) .SUMSQ(200) 

C  COMMON  STATEMENTS 

COMMON  IT1 . IT2.N0P.EN.ENBL0K, I  NT .M , NBLOK . NPT 1 .STBLOK . W I DTH . I  X  I N . EM 
1  .ANUM.CONTPL.NCAL.CONLOG 

COMMON  XMULT.XADD.YMULT.YADD  , AMPLOW , DBLOW . AMPH I , DBH I 
COMMON  IX.  I  TEMP.  AMP.  AVER. I CT . SD . SUMAMP . S ANUM . DB . A . SUMSQ 
C  PERFORM  LOOP  FOR  EACH  STATION 

DO  1  J=NPT 1 .M , I  NT 

C  IF  ABOVE  CURVE  DO  NOT  USE 

IF ( AMP( J )-77777. ) 10.1,1 

CHANGE  CROSS  SECTION  BACK  TO  DB/SQUARE  METER 
FOR  DISTRIBUTION 
10  AMP(J)=10.*LOG10F (AMP(J) ) 

ADD  CONSTANT  TO  GET  POSITIVE  NUMBER 

12  AMP(J)=AMP(J)+CONTPL 

IF  BELOW  MINIMUM  BOX  EDGE  COUNT  AS  SPECIAL  CASE 
IF(AMP(J)-STBL0K)2.3,4 

2  I  CT  < J.101 )  =  I CT ( J.101 )+l 
GO  TO  1 

3  ICT ( J. 1 ) =ICT ( J,1 ) +1 
GO  TO  1 

IF  ABOVE  MAXIMUM  BLOCK  EDGE  COUNT  AS  SPECIAL  CASE 

4  IF ( AMP ( J ) -E NBLOK ) 5 . 6 . 7 
7  ICT (J.102 )  =  I C  T (J,102)+l 

GO  TO  1 

6  ICT ( J, NBLOK) =ICT ( J.NBLOKJ+l 
GO  TO  1 

OTHERWISE  LOCATE  IN  PROPER  BOX  AND  ADD  ONE  TO  CURRENT 
VALUE  OF  THAT  BOX 

5  ANOB=AMP ( J ) /WIDTH+1 . 

14  NOB= ANOB 

13  ICT ( J , NOB )=ICT (J»N0B)+1 
1  CONTINUE 

RETURN 
END 
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n  n 


*  LIST 

*  LABEL 

*  SYMBOL  TABLE 

SUBROUTINE  STAT 

STAT  STATISTICAL  ANALYSIS  FOR  AMPLITUDE  DATA 
DIMENSION  STATEMENTS 

DIMENSION  A (50)  .AMP (220) , AVER ( 2 00 ) » DB ( 50 ) .  I C T < 200 » 1 04 ) » I  TEMP (220) 
1.IXI220)  . SANUM ( 200 )  .SD1200 )  .SUMAMP ( 200 ) .SUMSQ (  200 ) 

C  COMMON  STATEMENTS 

COMMON  I  T  1  .  I  T2  , NOP  .  EN  ,  ENBLOK  ,  I  N T  ,M  ,  NBLOIC  ,  NP T  1  , S TBLOK  ,  W  I  D TH  ,  I  X  I  N  ,  EM 
1  .ANUM.CONTPL.NCAL.CONLOG 

COMMON  XMULT  .XADD.YMULT  »YADD  ,  AMPLOW . DBLOW . AMPH I  , DBH I 
COMMON  IX  .  I  TEMP  .AMP  .AVER  .  I C T . SD . SUMAMP , SANUM , DB , A . SUMSQ 
C  FORMAT  STATEMENTS 

100  FORMAT (21H  ACCUMULATOR  OVERFLOW) 

101  FORMAT (3H0J=I5.18H  POINTS  TOO  LOW=I5) 

102  FORMAT (3H0J=I5.19H  POINTS  TOO  HIGH=I5) 

C  PERFORM  LOOP  FOR  EACH  STATION 

DO  50  J  =  NPT 1 .M , I  NT 

C  SET  SUMS  EQUAL  TO  ZERO  AT  START 

AVER ( J ) =0. 

SD ( J  )  =0  • 

c  sanum=number  of  points  to  use  in  analysis 

C  POINTS  ABOVE  CURVE  HAVE  BEEN  ELIMINATED 

SANUM ( J ) =ANUM-FLOATF (  I CT ( J  .  103 )  ) 

50  CONTINUE 

C  IF  ANY  OVERFLOW  HAS  OCCURRED  IN  ADDING— EXIT  HERE 

1  IF  ACCUMULATOR  OVERFLOW  3.4 

3  WRITE  OUTPUT  TAPE  3,100 
CALL  EXIT 

4  DO  5  J  =  NPT 1 »M , I  NT 

C  COMPUTE  AVERAGE  EACH  STATION 

AVER(J)=SUMAMP(J) /SANUM! J) 

IF (  ICT (J  .101 )) 17,18,17 

17  WRITE  OUTPUT  TAPE  3 , 1 0 1 , J , I CT ( J , 1 0 1 ) 

18  I F ( ICT (J.102 )  )  19,5.19 

19  WRITE  OUTPUT  TAPE  3 . 1 0 2  .  J  .  I C T ( J , 1 0 2 ) 

5  CONTINUE 

C  COMPUTE  STANDARD  DEVIATION 

6  DO  14  J  =  NPT 1 .M , I  NT 

SD ( J ) =SQRTF ( SUMSQ (J ) /SANUM ( J ) -AVtR ( J ) *AVER ( J ) ) 

14  CONTINUE 
RETURN 
END 
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* 

FAP 

COUNT 

51 

ENTRY 

READA 

ENTRY 

UNPAK 

ENTRY 

SKIP 

* 

STORE  INDEX  REGISTER  1 

SKIP 

SXD 

TEMPI, 1 

* 

PLACE  NUMBER  OF  RECORDS  TO  SKIP  INTO  ACCUMULATOR 

* 

and  into  index  REGISTER  1 

CL  A* 

1  *  4 

POX 

*1 

* 

GO  BACK  TO  MAIN  PROGRAM  IF  NSKIP=0 

TZE 

RETURN 

TM1 

LOOP3 

* 

SKIP  NSKIP  RECORDS  IN  FORWARD  DIRECTION 

IF  NSKIP 

L00P2 

TCOA 

* 

RTBA 

8 

T  I  X 

LOOP2  *1  *  1 

RETURN 

LXD 

TEMPI *1 

TRA 

2*4 

* 

bACSPACE  NSKIP  RECORDS  IF  NSKIP  - 

L00P3 

TCOA 

* 

BSRA 

8 

T  I  X 

LOOP3  *  1  *  1 

TRA 

RETURN 

* 

READ  BINARY  RECORD  FROM  A8  INTO  CORE 

READA 

RTbA 

8 

RCHA 

Cl 

TCOA 

* 

* 

UNPACK  TIME  CODE  WORDS  FROM  WORD  1  INTO 

I  T 1  AND 

LDQ 

=  0 

CLA 

I  T  1 

LGR 

18 

ALS 

18 

STO 

IT1 

STQ 

I  T  2 

* 

UNPACK  NUMBER  OF  POINTS  ON  TRACE  FROM  WORD  2 

CLA 

K 

LRS 

18 

ALS 

18 

STO 

K 

* 

RETURN  TO  MAIN  PROGRAM 

RET 

TRA 

1*4 

* 

PLACE  NEXT  WORD  TO  UNPACK  INTO  ACCUMULATOR 

UNPAK 

CLA* 

1*4 

LGR 

18 

ALS 

18 

* 

STORE  FIRST  18  BITS  INTO  IXl 

STO* 

2*4 

* 

STORE  LAST  18  BITS  INTO  1X2 

SLQ* 

3*4 

* 

RETURN  TO  CALLING  PROGRAM 

TRA 

4*4 

* 

STORAGE  FOR  RECORD  READ  IN  FROM  A8 

Cl 

IOCP 

I  T  1  *  *  1 

IOCP 

K*  *  1 

IOCPN 

C*  *2 
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I  ORP 

I  OCD 

c 

BSS 

TEMPI 

BSS 

TEMP2 

BSS 

I  T  1 

COMMON 

I  T  2 

COMMON 

K 

COMMON 

AAAA 

COMMON 

IX 

COMMON 

END 

IX-219* ,220 
0,0*0 
1 
1 
1 
1 
1 
1 

22 

220 
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APPENDIX  B 


SJOB 

SEXECUTE 
S I B JOB 

SIBFTC  PLOAC 


W I LSON ♦ FLORENCE 
I  B  JOB 
GCUSOURCE 
LIST  *DECK  *SDD 


CA1/3 

C 


21 

22 

23 

115 

111 


C 

C 

C 

C 

c 

c 


AMPLITUDE  PLOTS 


1/3/66 


1/7/66 


70 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


4020  PLOT  FOR  AMPLITUDE  PROGRAM 
DIMENSION  STATEMENTS 

DIMENSION  TlTLElt 12) *TITLE2(7)*TITLE3(10) *TITLE4(12) ♦ X 1 ( 100) *X4( 10 
10)  *  ST A (200)  9  AVDB ( 200 )  *AVDBP(200)  *AVDBM(200> 

FORMAT  STATEMENTS 
FORMAT ( 1415) 

FORMAT ( 12A6) 

FORMAT ( 14F5.0 ) 

FORMAT ( IX  » I5»24Xf3F12.5 ) 

FORMAT ( 1H, 12F10.3 ) 

STORE  LABEL  FOR  ID  FRAME 
CALL  STOIDVt 15HAMPLITUDE  PLOT  S  9  3  > 

SELECT  CAMERAS 
CALL  C AMRA V ( 2 ) 

MM=NUMBER  OF  TIME  INTERVALS  TO  PLOT 
NBLOK=NUMBER  OF  BLOCKS  IN  FREQUENCY  DISTRIBUTION 
M=NUMBER  OF  STATIONS  PER  TIME  INTERVAL 

IFFREQ=-1  IF  DISTRIBUTION  IS  ON  TAPE  BUT  NOT  TO  BE  PLOTTtD 
I FFREQ  =  0  IF  DISTRIBUTION  TO  BE  PLOTTED 

I F FREQ=+1  IF  DISTRIBUTION  IS  NOT  WRITTEN  ON  INPUT  TAPE 
READ(  5  9  21  )MM*NBLOK9M9 IFFREQ 
WRITE  (6921 )MM9NBLOK 9 M9 IFFREQ 

XM  I NFR  =  M I N I MUM  LIMIT  OF  DISTRIBUTION  BOXES 
LIMIT  OF  DISTRIBUTION  BOXES 
FOR  TIC  MARKS  FOR  X  AXIS  OF  DISTRIBUTION 
FREQUENCY  OF  PLOT 

FOR  TIC  MARKS  FOR  Y  AXIS  OF  DISTRIBUTION 
STATION  NUMBER  FOR  AVERAGE  PLOT 
STATION  NUMBER  FOR  AVERAGE  PLOT 
FOR  TIC  MARKS  FOR  X  AXIS  OF  AVERAGE 

Y  SCALE  VALUE  FOR  AVERAGE  PLOT 

Y  SCALE  VALUE  FOR  AVERAGE  PLOT 
FOR  TIC  MARKS  FOR  Y  AXIS  OF  AVERAGE 

READ  ( 5  *  2  3 ) XM I NFR 9XMAXFR 9 DXFR 9 YMAXFR 9 DYFR * XMI NA * XMAXA  9 DXA 9 YM I NA 9 
1 YMAX A  9  DY  A 

WRITE( 6*23)XMINFR9XMAXFR9DXFR,YMAXFR9DYFR9XMINA9XMAXA 9  DXA  9  YM I NA  9 
IYMAXA9DYA 
DO  100  LM  =  1 9 MM 
I  F ( IFFREQ ) 2039201 ,200 

IF  DISTRIBUTION  IS  ON  TAPE  BUT  NOT  TO  BE  PLOTTED 
SKIP  OVER  THIS  SECTION  OF  INPUT  TAPE 
20  3  READ (17 *22) ( T I TLE 1 ( J )  9 J- 1  ♦  12  ) 

READ ( 17*22)  < T I TLE2 < J >  ♦ J  =  1  *  7 ) 

READ ( 17  922  )  ( T I TLE 3 ( J )  9 J= 1  *  10  ) 

DO  204  J  =  1 9 M 

READ ( 17*22)(TITLE4(JJ)9JJ=l9l2) 

READ ( 17  *  111 )  ( XI ( JJ ) 9X4 ( JJ ) * JJ  =  1 9NBLOK ) 

204  CONTINUE 

CO  PLOT  AVERAGES 
GO  TO  200 


XMAXFR=MAXIMUM 
DXFR= I NCREMEN  T 

ymaxfr=maximum 
DYFR= INCREMENT 
XM  I  NA  =  M I N I MUM 
XMAXA=MAXI MUM 
DXA=INCREMENT 
YM I NA=M I N I  MUM 
YMAX A=MAX I  MUM 
DYA= INCREMENT 


READ  TITLES 
201  READ  (  17,22) (TlTLElt J) 
WRITE  (6,22) (TlTLElt J) 


FOR  DISTRIBUTION 
J=1 *12 ) 

J= 1  *  1 2  ) 


FROM  INPUT  TAPE 


c 

c 


c 

c 


c 

c 

c 

c 


c 


c 


c 


c 

c 

c 

c 

c 


c 


352 


401 

356 

357 


355 


READ  (  17*22)  (TITLE2U)  *J  =  1*7) 

WRITE  (6*22)  ( T I TLE2 ( J  ) *J=1*7) 

READ  ( 17*22) (TITLE3(J) *J=1*10) 

WRITE  (6*22) (TITLE3(J)*J= 1*10) 

ADVANCE  FRAME 
CALL  FRAMEV 

PRINT  OUT  TITLE  FRAME  FOR  TIME  INTERVAL 
CALL  PRINTV(72*TITLE1*100*500) 

CALL  PRINTV(39*TITLE2*100*450) 

L  =  0  FOR  FIRST  PLOT  THIS  FRAME 

L  =  C 

READ  DATA  FOR  NEXT  STATION  TO  PLOT 
DO  354  J  =  1  *  M 

READ  ( 17*22) (TITLE4(JJ)*JJ=1 *12) 

WRITE  (6*22)  ( T I T  LE4 ( J J ) ,JJ=1,12) 

READ  ( 17*111)  (XI (JJ) *X4(JJ>  ,JJ=l*NBLOK) 

IS  THIS  FIRST  PLOT  THIS  FRAME 
IF (L)355*356*355 

YES - SET  CODE  FOR  NEXT  STATION 

L=  1 

ADVANCE  FRAME 
CALL  FRAMEV 

PRINT  OUT  IDENTIFYING  TITLES  AT  TOP  OF  FRAME 
CALL  PRINTV(72*TITLE1*100*1010) 

CALL  PRINTV ( 39*T I TLE2*100*994) 

CALL  PRINTV(40,TITLE3*100*978) 

CALL  PRINTV(68*TITLE4*100*950) 

GET  SCALE  FACTORS  FOR  TOP  PLOT 
CALL  XSCALV( XMINFR  * XM AXFR ♦ 1 00  *  50 ) 

CALL  YSCALV (0.,YMAXFR *525*98) 

OUTLINE  TOP  PLOT 
CALL  LINEV( 100*525*100*925) 

CALL  LINEV( 100,525*973*525) 

CALL  LINEV(973,525*973,925) 

CALL  LINEV( 100,925,973,925) 

LABEL  AXES  TOP  PLOT 

CALL  LINRV( 1 *500*520*525  * XMINFR  *XMAXFR , DXFR *  1  *  1 ♦ 5  * 8 ) 

CALL  LINRVt 2 *50*95 *100*0. ♦ YMAXFR *DYFR , 1 ♦ 1 *4  *  10 ) 

CALL  PR  I  NT V (  36*36HCROSS  SECTION  (DB  ABOVE  1  SQ*  ME T ER ) , 390  * 480 ) 
CALL  APRNT V ( 0  *- 14  *  15,15HNUMBER  OF  CASES  * 2S , 850 ) 

SET  LOWER  Y  LIMIT  OF  PLOT 

I Y2  =  52  5 

GO  TO  PLOTTING  ROUTINE 
GO  TO  360 

NO - SET  CODE  FOR  ADVANCING  FRAME  NEXT  TIME 

L=0 

GET  SCALE  FACTORS  FOR  LOWER  PLOT 
CALL  YSCALV(0. * YMAXFR ♦ 35  *  588  ) 

OUTLINE  LOWER  PLOT 
CALL  LINEV( 100*35,973,35) 

CALL  LINEV( 100*35*100*435) 

CALL  LINEV( 100*435*973*435) 

CALL  LINEV(973,35,973,435) 

LABEL  AXES  LOWER  PLOT 

CALL  L I  NR V ( 1  *  460  * 435  *  440  * XM I NFR , XMAXFR *DXF R , 1  *  1  *  5  * 8 ) 

CALL  L I NRV(2* 50 *95* 100*0. , YMAXFR *DYFR, 1,1 *4,10) 

CALL  APRNT V (0,-14,  15*15HNUMBER  OF  CASES , 25 *360 ) 

CALL  PRINTV(68*TITLE4*100,20) 


52 


n  o 


c 

c 


c 


c 


c 


c 


c 


c 


c 


SET  LOWER  Y  LIMIT  OF  PLOT 

I  Y2  =  35 

PLOT  FREQUENCY  DISTRIBUTION  FOR  EACH  BLOCK  THIS  STATION 
360  DO  367  II  =  1,NBL0K. 

IFIX4I  II ) >367,367,340 
340  AA A  =  X 1 (II) 

IX1=NXV ( AAA ) 

AAB=X4 (II) 

I Y1=NYV ( AAB ) 

CALL  L  INEVI  1X1 , IY1 , 1X1 , IY2 ) 

IX1=IX1+1 

CALL  LINEVI  1X1  ,IY1  ,1X1  ,IY2 ) 

1X1=  1X1-2 

CALL  LINEVI  I  X 1  ,  I Y 1  ,  I  X 1  ,  I Y2 ) 

367  CONTINUE 

GO  PLOT  DISTRIBUTION  OF  NEXT  STATION  IF  ANY 

354  CONTINUE 

BEGIN  AVERAGE  PLOT 
ADVANCE  FRAME 
200  CALL  FRAMEV 

READ  TITLES  FROM  INPUT  TAPE 
READ  (17,22)  (TITLEl(JJ)  ,JJ=1, 12) 

WRITE(6.22)(TITLE1(JJ),JJ=1,12) 

READ  (17,22)  (TITLE21JJ)  ,JJ=1,7) 

WRITE (6,22) (TITLE2( JJ> »JJ=1»7) 

PRINT  TITLES  FOR  AVERAGE  PLOT 
CALL  PRINTVl 72 *T ITLE1 *  100, 1000) 

CALL  PRINTV(42,TITLE2,100,975) 

CALL  PRINTVl  28.28H  AVERAGE  PULSE  CROSS  SEC T I  ON  ,  100 , 95  0 ) 

SET  SCALE  FACTORS  FOR  AVERAGE  PLOT 
CALL  XSCALV(XMINA,XMAXA,100,50) 

CALL  YSCALV(YMINA,YMAXA, 100,123) 

OUTLINE  PLOTTING  AREA 
CALL  LINEVI 100,100,973,100) 

CALL  LINEVI 100,100,100,900) 

CALL  LINEVI 100,900,973,900) 

CALL  LINEV1973, 100*973, 900) 

LABEL  AXES 

CALL  LINRVI 1,80,95,100,XMINA,XMAXA,DXA,1 ,1,4,8) 

CALL  LINRV(2,50,95,100,YMINA,YMAXA,DYA,1 ,1 ,5,10) 

CALL  PRINTVl  1 1  , 1 1 HST AT  I  ON  NO. ,450, 25) 

CALL  APRNT V (0,-14,  34.34HCR0SS  SECTION  DB  ABOVE  1  SQ.  METER, 20, 780 
1  ) 

PLOT  AVERAGES  AND  STANDARD  DEVIATIONS  EACH  STATION 
DO  135  J J= 1 , M 

READ  I  17,115) JK.AVDBI JJ) ,AVDBP( JJ) , AVDBM I J J ) 

135  ST A I J J ) = JK 

DO  361  J  =  1 , M 
XAA  =  ST  A  I J ) 

Y  AA  =  AVDB I J ) 

CALL  P0INTV(XAA*YAA,1 ) 

AAA= AVDBP I J ) 

I Yl  =  NYV I  AAA ) 

I  X 1  =  NX V  I  X AA ) 

1X2= 1X1+5 
IX3=IXl-5 

CALL  LINEVI  1X3, IY1, 1X2,  I Y 1  ) 

AAA= AvDBM I J ) 
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IY2=NYV(AAA) 

CALL  LINEVf  1X3, IY2. 1X2, IY2) 
I  Y4  =  NY V ( Y  AA ) 

I Y  5= I Y  4+5 
I Y6= I Y  4-5 

CALL  LINEVf IX1.IY5.IX1.IY1) 
CALL  LINEVf IX1.IY2.IX1.IY6) 


c 

GO  DO  NEXT  STATION 

■ 

361 

CONTINUE 

c 

PLOT  NEXT  TIME  INTERVAL  IF  ANY 

100 

CONTINUE 

GO  TO  70 

M 

END 

*  SUBROUTINE  TO  DEFINE  INPUT  TAPE 


S1BMAP 

UNI  7 

6 

ENTRY 

•  UN  17* 

•  UN 1 7 • 

PZE 

UNI T17 

UN  I T 1 7 

FILE 

END 

*B<2 ) *HOLD* BLOCK- 22 
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